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FOREWORD 


This  report  was  prepared  by  Analytical  Services  &  Materials,  Inc.,  Hampton  Virginia  for 
WL/FIBEC,  Wright-Patterson  Air  Force  Base,  Ohio  under  contract  F33615-94-D-3212, 
“Structural  Integrity  Analysis  and  Verification  for  Aircraft  Structures.”  The  U.S.  Air 
Force  project  engineer  was  Mr.  James  A.  Harter,  WL/FIBEC.  The  government  contract 
manager  was  Lt.  Davis  S.  Conley.  The  period  of  performance  for  this  report  was 
November  1995  through  July  1996. 

The  work  performed  under  report  (Delivery  Order  0006)  was  performed  by  Analytical 
Services  &  Materials,  Inc.  personnel  located  at  the  WL/FIBE  Fatigue  &  Fracture  Test 
Facility,  Bldg.  65,  Area  B,  Wright-Patterson  AFB,  OH  and  their  subcontractors.  The 
Principal  Investigator  of  this  research  was  Mr.  Kevin  L.  Boyd.  The  authors  of  this  report 
were  Mr.  Kevin  L.  Boyd,  Mr.  Srinivas  Krishnan,  Mr.  Alex  Litvinov,  Dr.  Mohan  Ratwani 
(R-Tec;  Sections  2.0  &  3.0)  and  Dr.  Gregory  Glinka  (SaFFD;  Sections  4.0,  5.0  and 
Appendix  A).  Additional  technical  input  and  assistance  was  provided  by  Mr.  James  A. 
Harter. 


1. 


INTRODUCTION 


This  purpose  of  this  research  was  to  develop  advanced  structural  integrity  analysis 
methods  for  aging  aircraft  structures.  The  primary  focus  of  this  study  was  to  develop  a 
method  for  performing  fatigue  crack  growth  analyses  under  a  composite  (repair  patch)  on 
a  metallic  structure.  Also,  it  was  envisioned  that  these  methods  would  be  incorporated 
into  the  fatigue  crack  gro\vth  software,  AFGROW,  so  that  this  promising  technology  can 
be  transitioned  to  U.S.  Air  Force  Logistic  Centers,  U.S.  Air  Force  contractors  and  other 
government  agencies  as  quickly  as  possible. 

Fatigue  cracks  growing  from  fastener  holes,  or  multiple  fastener  holes,  also  are  common 
problems  found  in  aging  aircraft.  To  address  these  problems,  analytical  methods  (weight 
functions)  were  developed  to  assist  engineers  in  the  analysis  of  fatigue  cracks  growing  in 
arbitrary  stress  fields,  commonly  found  around  cold-worked  fastener  holes  and  stop- 
drilled  cracks.  Furthermore,  an  approach  to  evaluate  multiple  cracking,  or  multiple-site 
damage  (MSD),  was  developed  as  a  future  enhancement  to  the  AFGROW  program. 

The  development  of  enhanced  damage  tolerance  methods  are  necessary  to  accurately 
determine  the  inspection  intervals  and  safety  of  flight  of  repaired,  aging  aircraft 
structures.  Current  methods  used  to  evaluate  the  structural  integrity  of  repaired  structures 
include  time  consuming  (and  limited)  finite  element  analyses  and/or  static  strength  design 
criteria.  While  these  methods  are  often  adequate  for  estimating  the  equivalent  stiffiiess 
and  state  of  stress  of  a  repaired  structure,  they  do  not  include  load  interaction  effects, 
cyclic  debonding,  or  give  an  adequate  estimation  of  the  remaining  life  (or  inspection 
intervals)  of  a  repaired  structure.  The  primary  questions  regarding  bonded  composite 
repair  patches  based  on  comments  from  ASIP  managers  at  USAF  Air  Logistics  Centers 
[1]  are:  (1)  is  the  composite  repair  viable?;  if  so,  (2)  how  long  can  I  keep  the  aircraft 
flying  with  this  repair?,  and  (3)  how  often  should  I  inspect  the  repair?  The  statements  are 
not  necessarily  made  to  underestimate  the  importance  of  finding  qualified  personnel  to 
perform  repairs  and  trained  personnel  with  proven  NDE  methods  to  evaluate  the 
structural  integrity  of  the  patch.  In  most  cases,  the  decision  to  repair  or  replace  structural 
components  must  be  made  quickly  so  that  the  aircraft  may  be  returned  to  service  as  soon 
as  possible.  This  requirement  precludes  the  use  of  time  consuming,  detailed  finite 
element  analyses  or  rigorous  testing  programs. 

Presently,  engineers  at  the  ALCs  do  not  have  integrated  analysis  tools  to  assist  in  the 
design  of  repairs  and  subsequently  perform  Damage  Tolerance  Analyses  (DTA)  to  assign 
inspection  intervals  to  repaired  structures  [1].  Therefore,  it  is  important  that  advanced 
analysis  methods  be  developed  and  implemented  to  assess  the  structural  integrity  of 
repairs  being  performed  at  various  US  Air  Force  Air  Logistic  Centers  (ALCs).  It  is  also 
important  that  analysis  methods  are  user  friendly,  accurate,  fast,  economical,  upgradable, 
well  documented,  and  be  supported  by  knowledgeable  USAF  personnel. 
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This  research  does  not  intend  to  discourage  the  use  of  finite  element  methods,  static 
strength  analysis  methods  or  rigorous  testing  programs  to  accurately  characterize  the 
structural  integrity  of  a  patched,  repaired  structure.  What  this  research  wants  to 
emphasize  is  that  composite-patched  repaired  structures  should  be  analyzed  with  methods 
capable  of  characterizing  the  dynamic  aspects  of  a  crack  growing  under  a  composite 
patch  in  a  repaired  structure.  The  USAF  Aircraft  Structural  Integrity  Program  (MIL- 
STD-1530A)  calls  for  aircraft  to  meet  the  damage  tolerance  criteria  found  in  MIL-A- 
83444  (Airplane  Damage  Tolerance  Requirements).  It  is  the  belief  of  the  author  that  the 
design  of  repairs  made  to  aging  aircraft  should  also  be  evaluated  under  damage  tolerance 
requirements  to  whatever  extent  possible. 

With  regard  to  composite  repaired  structures,  the  potential  of  repairing  cracked  metallic 
structures  with  composite  patches  for  cost  reduction  and  improved  structural  efficiency  is 
becoming  evident  from  the  studies  of  various  investigations  [2-9].  The  in-service 
application  of  the  repair  concept  are  found  in  the  repair  of  C- 141  weep  hole  cracks,  C-5 
fuselage  skin,  F-16  fuel  access  hole  in  lower  wing  skin  and  T-38  lower  wing  skin  at  ‘D’ 
panel.  In  the  majority  of  cases  where  damage  metallic  structures  have  been  repaired  with 
composite  patches,  conventional  mechanically-fastened  repairs  could  not  be  used  due  to 
either  non-availability  of  space  to  drill  holes  or  detrimental  effect  of  these  holes  on  the 
structural  integrity  of  the  repaired  structure.  The  composite  patch  repair  concept  offers  an 
excellent  potential  to  repair  structures  which  cannot  be  repaired  with  conventional 
methods  and  thereby  reduce  the  operational  cost  of  in-service  fleet. 

This  repair  process  essentially  involves  bonding  of  fiber-reinforced  composite  material 
over  the  cracked  metallic  material,  thereby  reducing  the  stress  intensity  of  the  crack.  This 
results  in  a  slowing  of  crack  growth  in  cracked  metal  parts  or  even  complete  retardation 
in  certain  cases  [4].  The  use  of  composite  patch  repair  has  several  advantages  over 
conventional  repair,  namely: 

1.  More  efficient  load  transfer  from  the  cracked  part  to  the  repair  patch,  resulting 
in  more  reduction  in  stress  intensity  factors  at  crack  tips  as  compared  to 
conventional  mechanically  fastened  repairs. 

2.  No  additional  stress  concentrations  due  to  drilling  of  holes  for  fasteners  as  in 
the  case  of  conventional  repairs. 

3.  Good  durability  under  cyclic  loading. 

4.  High  directional  stiffness  of  a  patch  allows  use  of  thin  patches. 

5.  Can  be  used  to  repair  curved  surfaces  and  complex  geometries. 


In  order  to  assure  the  flight  safety  of  structures  repaired  with  composite  patches,  it  is 
necessary  to  develop  analysis  methods  to  predict  the  life  of  repaired  structures.  Also,  the 
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analysis  techniques  should  be  easily  used  by  the  engineers  involved  in  designing  repairs 
■without  causing  any  delay  in  aircraft  do'wn  time. 

This  report  has  been  broken  down  into  three  sections;  development  of  composite  repair 
design  criteria,  development  of  damage  tolerance  methods  for  repaired  structures,  and 
development  of  weight  functions  and  approach  to  MSD.  The  development  of  composite 
repair  design  criteria  is  described  in  Section  2  and  development  of  damage  tolerance 
methods  for  composite  patch  repaired  structures  is  described  in  Section  3.  The 
development  of  weight  functions  is  described  in  Section  4  and  approach  to  MSD  is 
reported  in  Section  5. 
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2. 


REPAIR  DESIGN  CONSIDERATIONS 


A  number  of  factors  need  to  be  considered  in  designing  composite  patch  repairs  for 
metallic  structures.  From  the  point  of  a  structural  designer  the  following  factors  need  to 
be  considered: 

1 )  Repair  patch  material  and  thickness 

2)  Repair  patch  ply  orientation  and  stacking  sequence 

3)  Patch  configuration 

4)  Adhesive  material  and  thickness 

5)  Adhesive  disbonding 

6)  Effect  of  moisture 

7)  Metallic  material  thickness 

8)  Effect  of  thermal  stresses 

Consideration  should  be  given  to  the  location  of  the  repair  when  selecting  the  patch 
material,  configuration,  taper.  External  aircraft  surfaces,  particularly  in  areas  critical  to 
aerodynamic  smoothness  requirements  (such  as  wing  leading  edges)  require  coordination 
with  the  proper  authority  for  changes  to  external  aerodynamic  surfaces.  A  repair  to  a 
substructure  will  be  dictated  by  the  proximity  of  adjacent  structural  components.  Some 
of  the  important  factors  influencing  a  composite  patch  repair  will  be  discussed  in  the 
following  sections. 


2.1.  Repair  Patch  Material  &  Thickness 

A  number  of  composite  materials  are  available  for  bonded  repair  applications;  however, 
the  most  commonly  used  materials  are  boron,  graphite,  and  sometimes  GLARE.  Each  of 
these  materials  has  an  application  in  composite  patch  repair  depending  on  the  metallic 
material  thickness,  type  of  load  spectrum,  and  the  stress  level  in  the  spectrum.  For  lightly 
loaded  and  thin  structures  such  as  in  fuselage  structures  GLARE  may  be  suitable.  It  has 
the  advantage  that  it  is  relatively  cheap  compared  to  boron  and  graphite.  Also,  from  the 
design  point  of  view  GLARE  has  nearly  the  same  coefficient  of  expansion  of  the 
aluminum  base  material  and  produces  no  residual  compressive  stresses  during  the  in- 
service  operational  environment  of  the  aircraft.  This  can  be  an  advantage  in  certain  cases. 
However,  for  highly  loaded  structures  and  structures  experiencing  severe  spectrum 
loading  it  is  desirable  to  use  materials  such  as  boron  or  graphite.  Both  of  these  materials 
provide  greater  stiffness  when  compared  to  GLARE.  Boron  patches  are  generally  thinner 
than  graphite  because  boron  has  a  modulus  which  is  about  50%  higher  than  that  of 
graphite.  However,  the  cost  of  Boron  is  about  20  times  that  of  graphite. 

It  is  desirable  to  keep  patch  thickness  to  a  minimum,  as  thicker  patches  will  result  in  out- 
of-plane  bending  when  the  repair  patch  is  applied  to  one  side  only.  The  out-of-plane 
bending  causes  additional  tensile  stresses  in  the  metal  sheet  and  increases  the  stress 
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intensity  factors,  thereby  reducing  the  effectiveness  of  the  patch.  In  most  cases,  there  is 
an  optimum  patch  thickness  which  reduces  the  state  of  stress  at  the  damage  in  the  metal. 
Increasing  the  thickness  beyond  the  optimmn  can  result  in  significant  out-of-plane 
bending  and  reduce  the  effectiveness  of  the  patch. 

In  Reference  10,  guidelines  have  been  provided  for  repair  patch  thickness  for  two  sided 
and  one  sided  repairs  shown  in  Figure  1.  In  the  reference  the  repair  patch  thickness  tj  is 
recommended  as  follows: 

For  a  two  sided  repair  patch  (Figure  la),  thickness  tj  of  patch  on  each  side  is  given  by; 

(E,t,/2E,)<tj<(l.2E,t,/2E,) 

For  a  one  sided  repair  patch  (Figure  lb),  thickness  of  the  patch  is  given  by: 

(E,t,/E,)<t,<(l.2E,t,/E,) 

The  patch  thickness  is  intended  to  restore  original  stiffness  over  the  defect.  It  is  desirable 
to  provide  a  patch  thickness  to  restore  original  stiffiiess  of  the  structure  being  repaired. 
However,  a  stif&iess  of  up  to  20%  higher  than  that  of  base  metal  is  considered  acceptable. 
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Repair  Patch  Thickness  t2,  Modulus  in  Loading  Direction  Ey 
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Metal  Thickness  ti,  Modulus  Ei 

(a)  Cracked  Metal  Sheet  with  Repair  Patch  on  Both  Sides 


Repair  Patch  Thickness  t2,  Modulus  in  Loading  Direction  Ey 


^  Crack - ►[  ^  Adhesive 
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''Metal  Thickness  ti,  Modulus  Ej 

(b)  Cracked  Metal  Sheet  with  Repair  Patch  on  One  Side 


Figure  1.  Repair  Configurations 


2.2.  Repair  Patch  Ply  Orientations  and  Stacking  Sequence 

It  is  desirable  to  use  highly  orthotropic  patches  having  a  high  stiffness  in  the  direction 
normal  to  the  crack,  but  with  some  fibers  in  directions  at  45°  and  90°  to  the  primary 
direction.  This  configuration  has  been  found  to  prevent  matrix  cracking  under  bi-axial 
loading  and  in-plane  shear  loads  which  exist  for  typical  applications.  A  good  design 
practice  is  to  use  about  70%  0°,  plies  20%  45°/-45°  plies  and  10%  90°  plies.  A  typical 
repair  patch  is  shown  in  Figure  2.  When  graphite/epoxy  is  used  as  a  repair  material,  a 
barrier  ply  of  ‘S’  glass  is  used  between  aluminum  substrate  and  graphite  patch  to  prevent 
galvanic  corrosion.  A  wrap  of  181  woven  glass  is  sometimes  applied  on  the  patch  surface 
as  shown  in  Figure  2  to  provide  a  smooth  ply  transition  on  exterior  surface  repairs.  This 
has  additional  advantage  that  during  any  in-service  damage  woven  glass  is  damaged  and 
the  repair  may  still  be  intact. 


6 
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Figure  2.  Typical  Ply  Orientations 


The  stacking  sequence  in  a  repair  patch  is  also  an  important  consideration  in  repair  patch 
design.  It  is  desirable  to  put  +/-45  degree  plies  on  the  outside  surface  not  in  contact  with 
damage  area  and  zero  degree  plies  within  the  patch  laminate.  This  has  two  advantages,  1) 
during  any  in-service  impact  damage  +/-45  degree  plies  will  get  damaged  and  the 
structural  integrity  of  the  patch  is  not  significantly  reduced,  and  2)  the  presence  of +/-  45 
degree  plies  reduces  the  inter-laminar  stresses  at  the  edges  of  the  patch  and  thereby 
reduces  the  tendency  of  the  patch  to  delaminate  at  the  edges  during  in-service  cyclic 
environment.  For  thin  patches,  the  stacking  sequence  may  not  have  any  significant 
influence  on  load  transfer  to  the  patch.  While  this  may  often  be  true  for  thin  patches,  it 
may  not  apply  to  thick  patches,  therefore,  it  is  desirable  to  place  0  degree  plies  closer  to 
the  metal  surface  for  thick  repair  patches.  However,  consideration  must  also  be  given  to 
maintaining  balanced  laminate  symmetry  to  avoid  warpage  in  the  pre-cure  of  the  patch. 


2.3.  Repair  Patch  Configuration 

The  repair  patch  configuration  may  vary  and  is  primarily  dependent  on  the  structure’s 
repair  location  and  loading  direction.  Composite  patch  repairs  to  the  external  skin  are  a 
function  of  geometry  and  may  take  on  the  appearance  of  a  circular  or  elliptical  shape  with 
circumferentially  tapered  edges  (Figure  3),  and  are  dictated  by  the  aerodynamic 
smoothness  criteria.  Where  geometry  dictates,  the  prestaging  of  the  patch  is 
recommended. 
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Repairs  made  to  substructural  components,  such  as  spar  caps  and  rib  flanges  with 
cracking  at  fastener  hole  locations,  may  require  a  staged  or  precured  patch  as  shown  in 
Figure  4.  The  ends  of  the  patch  are  tapered  for  smooth  load  introduction  and  to  reduce 
the  peel  stresses  at  the  ends. 

To  efficiently  transition  the  load  into  the  repair  patch,  plies  should  be  stepped  at  0.1  inch 
for  every  0.010  inch  thick  stacking  of  plies.  This  taper  rate  applies  to  both  boron/epoxy 
and  graphite/epoxy  patches. 


Pt\'MA.P,V 

VoDA,0 


Figure  4.  Substructure  Repair 
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2.4.  Adhesive  Material  &  Thickness 


Room  temperature  cure  adhesives  are  not  considered  suitable  due  to  service  temperature 
requirements  of  180°F  in  major  aircraft  repair  applications  and  long-term  durability 
needed  for  repairs.  Room  temperature  cure  adhesives  are  paste  adhesives  and  generally 
do  not  result  in  uniform  bond  line  thickness  in  the  repair,  thus,  affecting  the  load  transfer 
to  composite  patch.  Hence,  high  temperature  cure  adhesives  are  preferred.  A  3  5 OT  cure 
film  adhesive  is  not  considered  desirable  as  the  curing  at  such  a  high  temperature  is  likely 
to  cause  undesirably  high  thermal  stresses.  Also,  an  aluminum  structure  exposed  to 
350°F  temperature  will  experience  a  degradation  in  mechanical  properties.  A  250°F  cure 
temperature  film  adhesive  is  considered  suitable  for  the  composite  patch  repairs  of 
aluminum  structure. 

Ductile  adhesives  such  as  FM-73  are  preferred  over  brittle  adhesives  such  as  FM-400  due 
to  the  greater  tendency  of  the  brittle  adhesives  to  disbond  around  the  damage  area  during 
flight  loading.  This  adhesive  disbonding  reduces  the  load  transfer  effectiveness  to  the 
patch. 

Adhesive  thickness  should  be  kept  as  small  as  possible  so  in  order  to  develop  full  bond 
strength.  A  thick  adhesive  bondline  causes  less  load  transfer  to  the  repair  patch  and 
increases  the  stress  intensity  factors  in  the  cracked  metallic  structure  and  crack  growth 
rate.  Successful  bonds  are  achieved  using  one  ply  of  0.01 -inch  thick  adhesive.  First,  the 
adhesive  is  allowed  to  thaw  to  room  temperature  after  removing  it  from  the  refngerator. 
Next,  the  adhesive  is  cut  to  the  size  of  the  patch  and  installed  on  the  base  material  bond 
surface  in  the  prescribed  location.  Lastly,  the  appropriate  procedure  for  applying  pressure 
and  following  the  recommended  cure  cycle  must  be  followed  to  achieve  the  desired  bond¬ 
line  thickness.  A  desirable  bond-line  thickness  ranges  from  0.004-inches  to  0.007-inches. 


2.5.  Adhesive  Disbonding 

Disbonding  in  the  adhesive  layer  can  occur  during  the  fabrication  of  the  repair  patch, 
during  the  in-service  spectrum  loading  as  the  damage  in  the  metallic  propagates  due  to 
fatigue,  or  due  to  in-service  impact  damage.  The  overall  effect  of  the  adhesive 
disbonding  is  to  reduce  the  load  transfer  to  the  repair  patch  and  reduce  the  effectiveness 
of  the  patch.  The  disbonding  in  the  adhesive  layer  due  to  in-service  fatigue  loading  is 
generally  elliptical  in  shape  for  through  the  thickness  cracks  and  the  leading  edge  of  the 
disbond  generally  coincides  with  the  leading  edge  of  the  crack.  The  ratio  of  major  to 
minor  axis  of  the  disbond  depends  on  the  type  of  adhesive  being  used.  The  effect  of 
disbonding,  caused  by  crack  propagation  due  to  in-service  loads,  should  be  taken  into 
consideration  during  the  design  process.  However,  if  the  disbonding  occurs  due  to  in- 
service  impact  damage  or  fabrication  process,  it’s  effect  on  structural  integrity  will 
depend  on  the  location  of  the  disbond  with  respect  to  the  crack  plane. 
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2.6.  Effect  of  Moisture 


Composite  materials  have  the  propensity  to  absorb  moisture  from  the  environment  during 
in-service  usage.  The  presence  of  moisture  in  the  repair  can  influence  the  damage  growth 
in  the  repaired  structure.  The  effect  of  moisture  must  be  taken  into  consideration  in 
designing  composite  patch  repair  and  life  prediction.  The  maximum  moisture  content  in 
composite  patch  under  typical  environment  is  generally  about  1%.  The  presence  of 
moisture  can  cause  increase  in  crack  grovvlh  rate  due  to  two  reasons;  1)  The  presence  of 
moisture  increases  crack  growth  rate  in  the  metal  compared  to  room  temperature 
environment,  and  2)  the  presence  of  moisture  affects  the  load  transfer  to  the  patch.  The 
repair  patch  design  should  be  verified  by  testing  with  repair  patches  that  are  moisture 
conditioned  to  represent  realistic  environments. 


2.7.  Effect  of  Metallic  Material  Thickness 

The  effectiveness  of  composite  patch  repair  concept  depends  on  the  thickness  of  the 
metallic  structure  being  repaired.  In  thin  metallic  structures  (thickness  less  than  0.12 
inch)  repaired  with  composite  patches  on  one  face  or  both  the  faces  the  crack  growth  in 
the  metal  during  subsequent  flight  loads  is  fairly  uniform.  However,  in  thick  metallic 
structures  (thickness  greater  than  0.12  inch)  repaired  with  composite  patch  on  one  side 
only,  the  crack  growth  is  not  uniform  through  the  thickness.  Fatigue  crack  growth  is 
retarded  on  the  repair  patch  side  as  compared  to  the  side  with  no  repair  patch.  The  non- 
uniform  crack  growth  is  due  to  the  out-of-plane  bending  and  the  load  transfer  to  the  patch 
on  the  repaired  side  of  the  structure. 

In  thick  metallic  structures  with  repair  patches  on  both  the  faces,  the  crack  growth  may 
not  be  uniform  through  the  thickness.  The  crack  gets  retarded  more  on  the  surface  than  in 
the  central  portion  of  the  specimen.  This  is  because  there  is  more  load  transfer  into  the 
patch  in  the  surface  region  as  compared  to  the  central  portion  of  the  plate.  In  such  cases 
the  crack  may  exhibit  tuimeling  behavior.  The  effect  of  the  nonuniform  crack  growth 
should  be  taken  into  consideration  when  designing  composite  patch  repairs  for  thick 
structures.  The  existing  analytical  techniques  do  not  account  for  this  effect  and  need  to 
be  modified  before  yielding  realistic  results. 
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2.8.  Effect  of  Thermal  Stresses 


The  operating  temperatures  of  typical  transport  and  fighter/trainer  aircraft  generally  do 
not  cause  significant  thermal  stresses  to  warrant  thermal  stress  analyses.  However,  the 
mismatch  in  thermal  expansion  properties  between  a  composite  patch  and  repaired 
metallic  structure  can  cause  residual  thermal  stresses,  which  may  accelerate  crack  grov^h 
in  aluminum  structure.  The  residual  stresses  are  greater  in  a  graphite  repair  than  a  boron 
repair  because  of  a  larger  thermal  mismatch  between  graphite  and  aluminum  (as 
compared  to  aluminum  and  boron).  These  thermal  stresses  are  produced;  (1)  during  the 
bonding  of  composite  repair  patches  to  aluminum  and  (2)  during  in-service  operation. 
The  residual  stresses  produced  in  aluminum  during  the  bonding  process  are  compressive 
and  hence  have  a  beneficial  effect.  This  effect  is  generally  neglected.  In  operating 
environments,  at  temperatures  above  room  temperature,  the  aluminum  is  subjected  to 
compressive  stresses  and  the  “closure”  effect  is  beneficial.  However,  in  operating 
temperatures  below  room  temperature,  the  aluminum  is  subjected  to  tensile  stresses, 
which  will  reduce  the  crack  closure  effect  and  accelerate  damage  growth.  Therefore,  if  an 
aircraft  is  operating  at  -65°F  for  long  durations,  damage  growth  will  be  accelerated  and 
should  be  taken  into  consideration. 

In  the  mathematical  formulation  of  the  problem,  AT  is  the  temperature  difference  that 
produces  thermal  stresses  as  outlined  in  Section  3.4  of  the  analysis.  The  temperature 
difference  can  be  defined  by: 


AT  =  T„,-T„p 

where: 

Tref=  “Effective”  temperature  of  repair  system 
Top=  Temperature  of  current  operating  environment 


T,ef  is  a  temperature  necessary  to  quantify  the  thermal  residual  stresses  inherent  in  the 
metal  due  to  the  composite  bonding  process.  The  choice  of  reference  temperature  is 
highly  subjective  (difficult  to  measure  experimentally)  and  must  be  selected  \vith  caution. 
The  estimation  of  T„f  is  often  debated  within  the  repair  community  and  is  often  used  as  a 
“correction  factor”  used  to  fit  test  results.  For  example,  one  analyst  may  recommend 
using  T,ef  =  250°F  (a  typical  cure  temperature  for  film  adhesives  used  in  composite 
repairs),  and  estimate  the  T,ef  at  room  temperature  to  be  approximately,  250-70°F=  180°F. 
A  second  analyst  may  estimate  that  the  actual  temperature  at  which  the  adhesive  stops 
flowing  and  “locks  in”  the  residual  stresses  to  be  approximately  120T,  yielding  a  T^^f  of 
120-70°F=  SOT.  At  present,  there  is  no  definitive  answer  to  this  question.  The 
estimation  of  T^^f,  is  discussed  more  detail  in  the  following  paragraphs. 

In  summary,  temperature  differences  due  to  either  the  composite  bonding  process  and/or 
in-service  loading  can  produce  tensile  or  compressive  stresses  depending  on  the  situation. 
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To  illustrate  the  thermal  effects  of  a  composite  patch  on  a  metallic  structure,  the  two 
cases  specified  are  outlined  below: 

CASE  1:  Effect  of  Thermal  Stresses  due  to  Composite  Repair  Bonding  Process 

Thermal  stresses  are  produced  during  the  bonding  process  when  the  coefficients  of 
expansion  of  two  materials  are  different.  If  boron  or  graphite  is  used  as  repair  materials, 
thermal  stresses  will  be  produced  during  bonding  process.  All  of  the  repairs  on  test 
specimens  have  shown  that  the  test  specimens  are  bent  with  repair  patch  on  the  outside 
and  metal  on  inside.  This  indicates  that  the  metal  is  under  compression  and  the  patch  is 
under  tension.  Therefore,  the  stresses  produced  in  metal  are  compressive  and  will  tend  to 
reduce  stress  intensity  factor  to  some  degree.  This  effect  is  generally  neglected. 


CASE  2:  Effect  of  Thermal  Stresses  due  to  In-Service  Loading 
This  case  considers  a  composite  patch  repair  that  has  already  been  performed  on  a 
structure  and  the  structure  is  subjected  to  in-service  temperature  changes.  If  a  structure  is 
subjected  to  a  operating  temperature  T^p  which  is  higher  than  Tref?  then  metallic  structure 
will  try  to  expand  more  than  the  composite  repair,  and  will  be  restrained  by  the  composite 
repair  as  shown  in  Figure  5.  This  will  cause  compressive  stresses  in  the  metallic  structure 
and  tensile  stresses  in  the  repair  patch  as  shown  in  Figure  5.  Thus,  the  stress  intensity 
factor  due  to  thermal  stresses  (caused  by  a  temperature  increase)  should  be  subtracted 
from  the  stress  intensity  factor  due  to  mechanical  loading. 
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Figure  5.  Repaired  Structure  Subjected  to  Temperature  Increase 


In  the  case  where  a  structure  is  subjected  to  temperatures  lower  than  the  reference 
temperature,  T^^f  (e.g.  an  aircraft  flying  at  -65°F),  the  metallic  material  will  tend  to 
contract  more  than  the  repair  patch  (Figure  6)  but  will  be  restrained  by  the  composite 
patch.  This  will  cause  tensile  stresses  in  metal  and  compressive  stresses  in  the  composite 
as  shown  in  Figure  6.  Thus,  the  stress  intensity  factors  due  to  temperature  decrease 
should  be  added  to  the  stress  intensity  factors  due  to  mechanical  loading. 
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Metal  Contraction 


Metallic  Material 
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Equilibrium  Position 


_ 

—  Tension  in  Metal 

- 

L - 

^Compression  in  Composite 

Composite  Contraction 

Figure  6.  Repaired  Structure  Subjected  to  Temperature  Decrease 


It  should  be  noted  that  the  increase  or  decrease  in  stress  intensity  factors,  due  to  thermal 
stresses,  may  not  occur  at  the  same  time  as  the  maximum  stress  in  the  spectrum.  For 
example,  in  a  fighter  aircraft,  the  maximum  stress  occurs  during  combat  and  the 
temperatoe  is  not  likely  to  be  -65®F.  When  this  temperature  occurs  the  maximum  stress 
in  the  spectrum  is  likely  to  be  corresponding  to  Ig  or  2g  load,  which  is  likely  to  be  much 
less  than  design  limit  load  or  maximum  spectrum  load.  Therefore,  the  superposition  of 
the  thermal  and  mechanical  stress  intensity  factors  must  be  carefully  considered  in  order 
to  obtain  accurate  stress  intensity  factors  used  in  life  prediction  analyses. 

In  conclusion,  these  idealized  models  are  not  always  representative  of  the  situation 
encountered  in  real  structures.  In  real  aircraft  structure,  the  deformations  producing 
thermal  stresses  are  local  to  the  vicinity  of  the  repair.  Therefore,  the  influence  of  the 
substructure’s  constraint  will  depend  on  the  proximity  of  the  substructure,  relative 
stiffhess  of  the  substructure,  type  of  fastening,  etc.  Within  the  present  scope  of  this 
program,  it  was  impossible  to  develop  and  verify  the  analytical  techmques  used  to 
determine  the  effect  of  these  constraints. 
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3. 


METHOD  TO  ANALYZE  CRACKED  METALLIC  STRUCTURES 
REPAIRED  WITH  COMPOSITE  PATCHES 


The  analysis  of  a  cracked  metallic  sheet  with  a  bonded  repair  patch  has  been  studied  by 
various  investigators  [11-21].  In  References  17-19,  finite  element  techniques  have  been 
developed  to  analyze  a  cracked  metallic  sheet  with  a  bonded  composite  repair  patch.  The 
finite  element  method  is  not  a  very  convenient  tool  for  use  by  engineers  who  are  not 
experts  in  modeling  and  requires  special  training.  Therefore,  there  is  a  need  to  develop 
simple  methods  to  analyze  cracked  metallic  structures  repaired  with  composite  patches. 
These  techniques  should  be  easy  to  use  and  not  require  expert  knowledge. 

In  the  following  sections,  a  mathematical  method  of  analysis  is  described.  The  details  of 
the  approach  are  discussed  in  the  following  paragraphs. 


3.1.  Formulation  of  the  Problem 

The  problem  of  an  adhesively  bonded  composite  repair  patch  to  a  cracked  metallic  sheet 
is  shown  in  Figure  7.  A  fatigue  crack  of  length  2a  is  assumed  in  the  metallic  sheet. 
During  fatigue  cycling,  a  debond  may  develop  in  the  adhesive  layer.  The  debond  is 
assumed  to  be  elliptical  in  shape,  with  major  axis  along  the  crack  direction  and  the 
leading  edge  of  the  debond  coinciding  with  the  crack  tip. 
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The  analysis  of  the  cracked  metallic  structure  repaired  with  a  composite  patch,  shown  in 

Figure  7,  is  based  on  the  following  assumptions: 

1.  Composite  repair  patch  is  considered  as  an  orthotropic  material  with  principal 
directions  along  the  crack  plane  and  at  right  angles. 

2.  The  thickness  of  the  metallic  sheet  and  composite  patch  are  small  compared  to  the  in¬ 
plane  dimensions  so  the  structure  may  be  considered  under  plane  stress. 

3.  The  variation  of  the  stress  through  the  thickness  of  the  metal  sheet  and  composite 
patch  is  neglected. 

4.  Metal  sheet,  composite  patch,  and  adhesive  layer  are  linearly  elastic  materials. 

5.  The  adhesive  thickness  is  small  compared  to  the  metal  and  patch  thickness  so  that  the 
adhesive  may  be  treated  as  a  shear  spring. 

6.  The  surface  shear  transmitted  through  the  adhesive  acts  as  a  body  force. 

7.  The  dimensions  of  the  plate  and  patch  are  sufficiently  large  so  that  infinite  plate 
solutions  for  Green’s  functions  may  be  used. 


For  the  purpose  of  analysis,  the  structure  in  Figure  7  can  be  represented  by  the  structure 
shown  in  Figure  8;  where  the  structure  in  Figure  8a  has  no  crack  and  the  structure  in 
Figure  8b  has  a  crack  with  loading  applied  to  the  crack  surface.  The  solution  to  the 
problem  of  Figure  7  is  given  by  the  perturbation  problem  of  Figure  8b,  as  the  structure  of 
Figure  8a  has  no  crack  and  hence  no  stress  intensity  factor. 


a  a 


a  a 

(a)  Bonded  Sheet  with  No  Crack  (b)  Perturbation  Problem 

Figure  8.  Superposition  Technique  for  Cracked  Sheet  with  Patch 


In  the  perturbation  problem  of  Figure  8b,  let  u„  v,  and  Uj,  V2  be  the  x,  and  y  components 
of  the  in-plane  displacements  in  metal  sheet  and  composite  patch,  respectively,.  The  shear 
stresses  and  Xy  are  the  components  of  shear  stresses  in  the  adhesive.  The  continuity 
conditions  in  the  displacements  of  sheet  and  patch  give: 

U1-U2  =(h,/p3)  T, 

V, -V2  =(h,/p3)  Xy  (1) 

where,  is  the  shear  modulus  of  adhesive  and  h^  is  adhesive  thickness. 

The  perturbation  problem  of  Figure  8b  may  be  considered  as  a  metallic  sheet  with  body 
forces  and  crack  surface  loading  as  shown  in  Figure  9a  and  composite  patch  with  body 
forces  as  shown  in  Figure  9b.  The  body  forces  acting  on  the  cracked  sheet  are  Xj  and  Y, 
at  Zo  (Xg,  yj  in  x  and  y  directions,  respectively.  The  body  forces  acting  in  the  composite 
patch  are  Xj,  Y2  at  (x^,  yj  acting  in  x  and  y  directions,  respectively. 
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h-*\  h- 


(a)  Cracked  Sheet  w/Body  Forces 
and  Crack  Surface  Loading 


K 


X2 


Disbonded  Region 


Composite  (Ex,  Ey,  Vx,  Gxy) 


(b)  Composite  Patch  w/Body  Forces 


Figure  9.  Loading  on  Metal  Sheet  and  Composite  Patch  in  Perturbation  Problem 


The  body  forces  between  the  sheet  and  composite  patch  are  transferred  through  the 
adhesive  shear  stresses.  The  relationship  between  the  body  forces  and  shear  stresses  is 
given  by: 


X,=-(t,/h,),  Y,=-(T,/h,) 

X,=-(t./h,)l  Y,=-(T,/h,)  (2) 

The  stresses  and  displacements  in  terms  of  complex  functions  tj)  and  v)/  are  given  by 
Reference  22: 


=  2[<t)(z)+  ‘t>(z)] 

(3a) 

CTy  -  a,  +  2fT^y  =  2[(z  -  z)f  (z)-  (l)(z)+  V|;(z)] 

(3b) 

2p(u  +  tv)  =  K  J<j)(z)-  |tt/(z)dz  -  (z  -  z)(t)(z) 

0  0 

(3c) 
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3.1.1.  Displacements  in  Orthotropic  Composite  Repair  Patches 

In  the  pertmbation  problem  in  Figure  9b,  the  body  forces  Xj  and  Yj  are  the  only  forces 
acting  on  the  patch.  The  stresses  and  displacements  in  the  patch  can  be  obtained  from  the 
solution  of  the  problem  of  concentrated  body  forces  acting  at  any  arbitrary  location  Zq. 
This  solution  is  used  as  Green’s  function  and  integrated  over  the  bonding  region  ‘D’  to 
obtain  displacements  and  stresses  due  to  forces  acting  in  the  entire  bonded  region.  The 
Kolosov-Muskhelishvili  functions  for  the  orthotropic  composite  under  concentrated 
loads,  X2  and  Y2  acting  Zg,  are  given  by  [23]: 

*1^2  (^l)“  (^11^2  Cl2^2)^®§(^l  ~^i) 

X!fe)  =  (C2,X,  +C,,Y,)log(z,  -4,)  (4) 

where, 

Zk=x  +  s„y,  §,,=Xo+s^yo,  (k  =  l,2) 
z  =  X  +  iy,  Zq  =  Xo  +  lYo  (5) 

C„  =(1/27t0[s2  +S2  +Si  +SlS2S2(v^Ey/Ejl/[(s,  -S^^S,  -Sl)(l-S2/s,)] 

C,2  =  (1/2;:/)  S2S2  +S1S2  +S1S2  +v^]/[(si  -S2)(s,  -Si)(l -S2/S,)] 

C2,  =(1/27:0  [s,  +  Si  +S2  +SiSiS2(v,,Ey/Ej]/[(s2  -s,)(s2  -S2)(l-Sl/S2)] 

C22  =(1/27:0  SiSl  +S2S1  +S2S1  +V^]/[(s2  -S,)(s2  -S2)(1-Si/s2)]  (6) 

In  Equation  (6)  s,  and  S2  are  roots  of 

s‘+(E./G.,-2'’xy +E./E,=0  (7) 

for  which  Im(Sk)  >  0. 

The  stresses  and  displacements  in  the  plate  may  be  expressed  in  terms  of  (t)2  and  X2- 
Assuming  that  the  body  forces  Xj  and  Yj  are  continuous  functions  of  Zo(Xo,yo)  defined  in 
bonded  region  ‘D’,  the  displacements  at  a  point  in  the  composite  may  be  written  as: 

U2  (x,y)  =  [K,  1  (x,y;  x^  ,yo  )X2  (xo,yo )  +  K,2  (x,y;  Xo  ,yo )  Y2  (x^  ,y  „ )]  dxody  0 

V2(x,y)=  JJ^[K2,(x,y;  Xo,yJX2(xo,yo)+K22(x,y;  Xo,yo)Y2(xo,yo)]dXodyo  (8) 
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where  kernels  Kjj  (i,j  =  1,2)  are  given  by  the  Green’s  functions  (4). 


3.1.2.  Displacements  in  Cracked  Plate 


The  displacements  in  the  cracked  plate  consist  of  two  parts,  one  due  to  uniformly  applied 
stress  on  the  crack  surface  and  the  other  due  to  the  body  forces  acting  on  the  cracked 
plate.  The  displacements  due  to  uniform  applied  stress  on  the  crack  surface  are  given  by: 


v„ 


(x,y)  =  cytK-l)Re 

(x,y)  =  a|K  +  l)lm 


-2ylm 

-2yRe 


z/(z' 

z/(z' 


-  a-  y^  ]+  (l  -  k)x}/4p  =  a  fj  (x,y) 
-a'  y^]+  (l  -  K)y  j/4p  =  af2  (x,y) 


(9a) 

(9b) 


where, 

V,  is  Poisson  Ratio  of  plate  material 
p  is  the  shear  modulus  of  plate  material 


k  =  (3-v,)/(1  +  Vi) 

The  displacements  u,2  and  v,2  in  the  cracked  plate  due  to  body  forces  are  given  in 
Reference  21 .  These  displacements  are  given  by: 

2p(u,2  +  v,2 )  =  S  {-K[log(z  -  Zo )+ log(z  -  zo )]+  k[0,  (z,Zo  )+ 0]  (z,zo  )]/2  -  [0i  (z,Zo )+  k^0,  (z,zo  )]/2 
+  (k — 1)[k02  (z)  ~  02  (z)]/  2 + (z — Zo  )95  (z,Zq  )/  (z—Zq^ 

+  s[(z-Zo)/(z-zo)-1+k03(z,Zo)-03(z,Zo)-04(z,zJ+k04(z,zoJ 
+  Rigid  Body  Displacement  (10) 


where. 
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S  =  (X,+Y,)/[27i(1  +  k)] 

S  =  (X,-Y,)/[27t(l  +  K)] 

6|(x.xj=logzzo-a^  +(z,’-aT(2’ 

1/2 

II 

<N 

CD 

logz  +  (z^ 

63  (z,  Zq  )  =  (zo  -  Zo  )/2(z  -  zo )]  1  +  (z^  -  a^  1 

(zo^  -a^J^ {(a^  -zo^i 

)  =  [(z  -  z)/2(z  -  zo )  l(z)-  l(zo )/ (z'  - 

-f] 

05  (z,  Zo )  =  04  (z,  zo )-  (z  -  z)  j(zo )/  2(z^  -  a^  ^ 

p] 

l(z)=  7t 
j(z)=  7t 

'(z=-a>f-z] 

(11) 

Assuming,  the  body  forces  Xj,  and  Y,  are  continuous  functions  of  (xo,yo)  defined  in  the 
bonded  region  ‘D’,  and  using  equations  10  and  11  as  the  Green’s  functions,  the 
displacements  in  the  cracked  plate  due  to  body  forces  may  be  written  as: 

u,2  (x,y)  =  [h„  (x,y;  Xoyo )  X,  (x^  y^ )+  h,^  (x,y;  Xgyo)  Y,  (xo  yo )]  dxodyg 

Vi2  (x,y)  =  [h2i  (x,y;  Xoy  o )  Xj  (xo,yo )+  K2  (x,y;  Xoyo )  Y,  (xo,yo )]  dXodyo  (12) 

The  total  displacement  in  the  cracked  plate  may  be  written  as: 

u,(x,y)  =  Ui, +U12,  Vi(x,y)=v„+v,2  (13) 


3.1.3.  Integral  Equations  for  x,  &  Xy 

The  \  and  Xy  are  the  adhesive  shear  stresses  in  the  x  and  y  directions,  respectively.  The 
relationship  between  the  body  forces  and  shear  stresses  is  given  by  Equation  2.  Using 
displacements  u„  Uj  and  v„  Vj  in  Equation  1  gives  the  following  system  of  integral 
equations  for  the  unknown  shear  stresses  and  Xy: 

(ha/HaK(x,y)  +  llk(x.y  ;  Xo  ,yo )  +  T,  (xo,yo )+  k,2  (x,y;  Xo ,y 0 ) + Xy  (xo,y 0 )]  dXodyo  =  a  f,  (x,y) 
(ha/fla)T^y(x.y)  +  ;xo,yo)+x^(xo_yo)+k22(x,y;  Xo,yo)+Xy(xo,yo)]dxodyo  =af2(x,y) 

(14) 
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The  functions  f,  and  £3  are  given  by  Equation  (9)  and. 


k,]  =  H,]/h,  +  Ki]/h2  k]2  -  H,2/h,  +  Kj2/h2 

k2i  =1^21 /h]  +K2i/h2  k22  +^22/^2 

Kernels  K^j  (i,j  =  1,2)  which  have  logarithmic  singularities  are  known  and  are  square 
integrable  in  region  ‘D’. 


3.1.4.  Stress  Intensity  Factors 

The  stress  intensity  factors  in  plane  stress  problems  may  be  expressed  in  terms  of 
Kolosov  Muskhelishvili  function  <t)(z)  as  follows  [24]: 

k,  -  /k2  =,  lirn,  [2(x  -  af  ^  [a^  (x,0) -  (x,0)] 

lim^  2[2(z  -  a)f  ^  (|)(z)  (1 6) 

In  the  present  problem,  the  shear  component  of  the  stress  intensity  factor  is  zero  due  to 
the  symmetry  in  loading  and  geometry  and  only  cleavage  component  kj  is  present.  The 
cleavage  component  is  obtained  by  adding  the  effects  of  crack  surface  loading  Po  and  the 
body  forces  X,  and  Y,.  The  stress  intensity  factor  kj  is  given  by: 

k,=a(ar+||,[h,(x  0  ,yo )  X,  (xo  ,yo )  +  h2  (xo  ,y  0 )  Y,  (xg  ,yo )]  dxody^  (17) 

where,  h,(Xoyo)  and  h2(xo,yo)  are  the  cleavage  components  of  the  stress  intensity  factor 
due  to  the  concentrated  body  forces  X,  and  Y,,  respectively.  The  stress  intensity  factors 
are  given  by: 


k,  -  jk2  =  (l/a)'^'  {-S  [a  +  l(zo )]  /(a  -  Zo )  +  X  S  [a  + 1  (zo )]  /(a  -  Zq  )  -  s(zo  -  z^ ) 

[(a  +  l(zo ))/ (a  -  Zo )^  +  J  (zo )/(a  -  Zq  )]  (18) 

h,(Xo,yo)  and  h2(Xo,yo)  are  the  coefficients  of  X,  and  Y^  respectively  in  real  part  of 
Equation  (18). 
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3.2.  Solution  of  Integral  Equations 

The  integral  equations  given  by  Equation  (14)  are  the  Fredholm  type  and  may  be  solved 
by  collocation.  The  region  ‘D’  is  divided  into  smaller  cells  and  unknown  shear  stresses 
and  Ty  are  assumed  to  be  constant  in  each  cell.  In  the  integral  equations  the  integral  is 
replaced  by  summation  and  using  collocation  a  2N  X  2N  system  of  equations  are 
obtained.  The  equations  are  given  by: 

N  N 

(ha  /l^a  )  (Xi  ,  y  j  )+  £  £  [k,  1  (Xi  ,yj ;  Xop  ,y oq  )  T,  (xqp  ,y Oq  )+  k,2  (xj  ,y J  ;  Xop  ,yoq  )  Xy  (xqp  ,y Oq  )]  A,p  Ayq 

p=l  q=l 

=  ^fi(xi,yj)  (19a) 

N  N 

(ha  /l^a )  -Ty  (xi ,  y  j )+  £  £  k,  (xi  ,y  j ;  Xop  ,yoq )  X,  (xpp  (x;  ,y j  ;  Xop  ,yoq )  Xy  (xop  ,yoq )]  A^p  Ayq 

p=l  q=l 

=  f^fi(xi»yj)  (19b) 

i  =  l, . N,  j  =  l, . N 

The  kernels  in  Equation  (19)  have  logarithmic  singularities,  hence  the  singular  part  of  the 
kernels  is  evaluated  separately  in  the  closed  form.  In  the  actual  integration  a  telescopic 
grid  is  used.  A  typical  grid  used  in  integration  is  shown  in  Figure  10,  where  only  a  quarter 
of  the  integration  region  is  shown  because  of  symmetry.  The  grid  size  is  kept  small 
above  the  crack  plane  for  a  distance  of  about  half-crack  length  ‘a’,  as  the  shear  stresses 
are  high  in  this  region.  The  grid  size  is  increased  as  the  distance  from  the  crack  increases. 
Theoretically  the  integration  region  is  infinite,  however,  in  the  perturbation  problem  the 
shear  stresses  in  the  adhesive  decrease  rapidly  away  from  the  crack  plane.  Hence,  the 
integration  region  is  kept  finite  as  shown  in  Figure  10.  The  integration  region  is 
determined  so  that  the  stress  intensity  factors  are  not  significantly  affected.  The 
convergence  of  the  solution  depends  on  the  crack  length  and  elastic  properties  of  an 
adhesive,  cracked  sheet,  and  composite  patch. 
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Figure  10.  Grid  Used  in  Numerical  Integration  (Due  to  Symmetry  Only  1/4  of  Panel  Shown) 


3.3.  Effect  of  Bending 

The  analysis  of  a  cracked  sheet  repaired  with  a  composite  patch  (in  Section  3.1)  is  for 
extension  loading  and  it  is  assumed  that  the  patch  or  sheet  has  no  bending  stiffness.  It  is 
applicable  to  cases  where  the  repair  patches  are  applied  to  both  the  surfaces  of  the 
cracked  sheet  as  shown  in  Figure  11a.  This  type  of  repair  is  not  always  feasible,  for  in 
the  majority  of  cases  the  access  for  repair  is  available  on  one  side  only,  and  therefore  the 
repair  patch  may  be  applied  on  one  surface  only  (Figure  lib).  The  presence  of  a  patch  on 
one  side  only  will  produce  out-of-plane  bending  at  the  crack  plane.  This  bending  will 
produce  tensile  stresses  in  the  sheet  and  compressive  stresses  in  the  repair  patch.  The 
tensile  stresses  due  to  bending  in  the  sheet  will  produce  stress  intensity  factors  in  addition 
to  those  produced  by  the  applied  membrane  stress,  a.  The  stress  intensity  factors  due  to 
bending  should  be  added  to  those  produced  by  membrane  stresses.  A  simple  method  of 
calculating  stress  intensity  factors  due  to  bending  [11-12],  is  discussed  in  the  following 
paragraphs  for  the  composite  repair  patch  geometry  shown  in  Figure  lib. 
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(a)  Cracked  Sheet  w/Repair  on  Both  Surfaces  (b)  Cracked  Sheet  w/Repair  on  One  Surface 

Figure  11,  Repair  Patch  Conflgurations 

In  a  cracked  sheet  with  a  central  crack  of  2a  and  far  field  applied  stress  of  or,  the  load 
released  due  to  the  presence  of  the  crack  at  the  crack  plane  is  (2ahi)G  (area  x  stress).  This 
load  is  entirely  taken  by  the  stress  singularities  produced  ahead  of  the  crack  tips.  In  a 
cracked  sheet  repaired  with  a  composite  patch  shown  in  Figure  lib,  the  same  load 
released  due  to  the  presence  of  a  crack  (2ah]C)  at  the  crack  plane  is  partly  taken  by  the 
stress  singularities  ahead  of  the  crack  tips,  as  in  the  case  of  a  single  cracked  sheet,  and  the 
remainder  of  the  load  is  transferred  to  the  repair  patch  through  the  adhesive.  The  load 
transferred  to  the  patch  produces  bending.  Let  the  load  transmitted  to  the  patch  be  2ah,CT, 
(area  x  stress)  where  cr,  is  equivalent  stress  transferred  to  the  patch. 

Assuming  the  thickness  of  the  adhesive  to  be  small  compared  to  plate  thickness,  the 
moment  produced  by  the  load  transfer  between  the  cracked  sheet  and  patch  is  given  by: 

M  =  2ah,a,[(h, +h2)/2]  (20) 

where  (h,  +  h2)/2  is  the  distance  between  the  centerline  of  the  metal  sheet  and  repair 
patch. 

or 

M  =  ah](h, +h2)a,  (21) 
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Assuming  that  the  bending  is  resisted  by  the  plate  width  equal  to  the  patch  width  W,  the 
maximum  bending  stress  in  the  layer  is  given  by: 


where  y^jg^  is  the  distance  of  the  outer  sheet  edge  from  neutral  axis  and  I  is  the  equivalent 
moment  of  inertia  of  the  section.  For  computing  the  moment  of  inertia  and  y^^^  the 
thickness  of  the  adhesive  may  be  neglected  and  the  section  of  the  patch  resisting  bending 
moment  (Figure  10a)  may  be  represented  by  that  shown  in  Figure  12b.  For  computing 
the  moment  of  inertia  and  y^3^  the  equivalent  section  is  shown  in  Figure  12c,  where  it  is 
assumed  that  Ey  is  greater  than  E,  and  the  section  is  uncracked. 


i. 
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(a)  Section  of  Cracked  Sheet  w/Patch 


(b)  Idealized  Section  for  Bending  Analysis 


Equivalent  Section  of  Composite  Patch  (Thickness  hj,  Ey  is  Modulus  in  Loading  Direction) 
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(Thickness  hi) 
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(c)  Equivalent  Section  for  Bending  Stress  Analysis 
Figure  12.  Equivalent  Section  for  Bending  Stress  Analysis 


Substituting  Equation  (22)  in  (21)  gives: 

Gb  =[ah,(hi+h2)a,y„„]/I 


(23) 


3.3.1.  Load  Transfer  Factor 

The  load  transfer  factor,  Lf,  is  defined  as  the  ratio  of  the  load  transferred  (to  the  patch)  to 
the  load  released  due  to  the  presence  of  a  crack. 

Lf  =(2ah,a,)/(2ah,G)=a,/a  (24) 
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Substituting  Equation  (24)  in  (23)  gives: 


Gb  =[ah,(h, +h2)LfCTy„,Jl 

ab/a  =  [ah,(h, +h3)Lfy„3,]/I  (25) 

=  BC,  Bending  Correction  Factor 

Therefore  the  effective  stress  used  in  computing  the  stress  intensity  factors  is  increased 
by  an  amount  Cb. 


3.3.2.  Load  Transfer  Factor  Computation 

In  a  cracked  structure  with  no  repair  patch  the  stress  intensity  factor  may  be  written  as: 


K,  =CT(7ia)’/^f(a/W)  (26a) 

where,  a  is  half  crack  length  and  f(a/W)  is  finite  width  correction. 

Kg  is  the  stress  intensity  factor  in  the  cracked  structure  with  repair  patch  and  applied 
stress  of  g.  The  equivalent  stress  g^  acting  on  a  sheet  with  no  repair  patch  and  giving 
stress  intensity  factor  of  Kg  is  given  by: 


Kg=Gg{7raff(a/W)  (26b) 

Using  Equation  (26),  Kg  may  be  expressed  as: 

K.=a,K,/cj  (26c) 

or 

<y,  =  K,c/K.  (27) 

The  stress  transferred  to  the  repair  patch  is  the  difference  between  the  remotely  applied 
stress  in  the  repaired  structure  and  the  uniform  in-plane  stress  that  produces  the  same 
stress  intensity  factor  in  a  single  sheet  as  in  the  repaired  structure.  The  stress  transferred, 
G„  may  be  written  as: 


G,  =G-Gg 


(28) 
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Using  Equation  (27), 


a,  =a-K,a/K, 

The  load  transfer  factor  is  given  as: 

L,=a,/a  =  l-K,/K,  (29) 

is  the  stress  intensity  factor  in  the  adhesively  bonded  structure  obtained  from 
mathematical  or  finite  element  analysis.  is  the  stress  intensity  factor  in  a  single  layer 
structure  having  the  same  crack  length  as  the  repaired  structure  and  can  be  obtained  from 
any  handbook  on  stress  intensity  factors.  Substituting  this  value  of  load  transfer  factor  in 
equation  (25)  gives  the  bending  correction  factor,  BC: 


BC  =  ab/a  =  [ah,(h,  +h2Xl-Ka/K3)y,3,]/I  (30) 

The  moment  of  inertia  I  and  y„,ax  are  obtained  for  the  equivalent  section  shown  in  Figure 
12c.  In  the  equivalent  section  shown  in  the  figure,  Ey  is  the  modulus  of  the  repair  patch 
in  the  loading  direction. 


3.4.  Effect  of  Thermal  Stresses 

In  a  cracked  metallic  structure  repaired  with  composite  patches,  the  coefficients  of  linear 
expansion  of  the  cracked  metallic  sheet  and  composite  repair  patch  are  different,  and  this 
difference  in  the  coefficients  of  expansion  gives  rise  to  thermal  stresses.  These  thermal 
stresses  may  be  due  to  the  residual  stresses  produced  during  the  curing  process  or  due  to 
the  in-service  environment  experienced  by  the  repair.  The  fabrication  of  the  panels  with 
graphite  or  boron  repairs  have  shown  that  residual  stresses  produced  by  fabrication 
processes  are  generally  compressive,  have  beneficial  effects  on  crack  growth,  and  are 
generally  neglected.  However,  the  residual  stresses  produced  by  in-service  change  of 
temperature  may  produce  tensile  or  compressive  stresses.  In  general,  an  increase  in  in 
service  temperature  compared  to  room  temperature  will  produce  compressive  stresses  and 
hence  a  reduction  in  stress  intensity  factors.  While  a  reduction  in  in-service  temperature 
will  produce  tensile  stresses  and  an  increase  in  stress  intensity  factors,  it  may  be  noted 
that  no  residual  stresses  are  produced  if  the  repair  material  is  GLARE  as  the  coefficients 
of  expansion  of  aluminum  and  GLARE  are  about  the  same. 

A  simplified  thermal  stress  analysis  of  a  metallic  structure  bonded  to  composite  was 
presented  in  Reference  25.  This  stress  analysis  was  used  in  the  present  work  to  compute 
stress  intensity  factors  due  to  thermal  stresses.  Consider  the  bonded  structure  shown  in 
Figure  lib.  Let, 
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tti  =  Coefficient  of  expansion  of  metallic  material 

tty  =  Coefficient  of  expansion  of  composite  repair  patch  in  loading  direction 
AT  =  Temperature  difference  producing  thermal  stresses  (T„rTop) 

L  =  Dimension  of  the  patch  in  loading  direction  at  right  angles  to  crack  plane 

E,  =  Modulus  of  elasticity  of  metallic  material 

Ey  =  Modulus  of  elasticity  of  composite  material  in  loading  direction 

The  thermal  stresses  are  computed  as  follows: 


Expansion  of  metallic  material  =  a,T(L/2) 

Expansion  of  repair  patch  =ayT(L/2) 

Differential  deformation  producing  thermal  stresses  =  a]T(L/2)-  a  ,T{L/2) 

=  (a,-a,)AT(L/2)  (31) 

Strain  produced  by  thermal  stresses  =  (a,  -  )  AT(L/2)  =  (a,  -  a  )aT  (32) 

Residual  stress  in  metallic  sheet  due  to  thermal  mismatch  is  given  by, 

E,  (a, -tty) AT  (33) 

This  far  field  stress  is  used  to  compute  stress  intensity  factors  in  the  structure  with  a 
repair  patch.  It  may  be  noted  that  the  thermal  stress  analysis,  discussed  here,  is 
approximate  and  does  not  account  for  two  factors  namely;  1)  effect  of  adhesive 
deformation  in  relieving  strain  and  reducing  the  effect  of  thermal  mismatch,  and  2)  the 
effect  of  relative  differential  expansion  between  composite  and  metal  at  right  angles  to 
the  loading  direction  along  the  crack  plane. 


3.5.  Effect  of  Disbonding  in  Adhesive  Layer 

Experimental  results  [11-12]  have  shown  that  as  the  crack  propagates  under  fatigue 
loading,  disbonding  occurs  in  the  adhesive.  This  disbonding  carmot  be  predicted  by  a 
static  analysis  because  it  occurs  due  to  cyclic  loading  at  much  lower  stresses  than  those 
predicted  by  static  analysis.  The  prediction  of  disbonding  due  to  cyclic  loading  is  rather 
complex  and  was  not  attempted  in  the  present  research.  The  disbonding  in  the  adhesive 
layer  influences  the  load  transfer  to  the  repair  patch,  and  ultimately,  the  stress  intensity 
factors  in  the  cracked  sheet.  The  presence  of  a  disbond  will  cause  less  load  transfer  from 
the  cracked  sheet  to  the  repair  patch  and  thereby  increase  the  stress  intensity  factors  and 
hence  crack  growth  rate.  The  presence  of  cyclic  disbonding  should  be  taken  into 
consideration  in  making  any  life  predictions. 
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The  experimental  studies  were  used  to  account  for  the  disbond  in  the  adhesive  layer. 
These  studies  have  shown  that  the  disbonding  in  the  adhesive  layer  tends  to  be  elliptical 
in  shape  with  the  leading  edge  of  the  disbond  coinciding  with  the  leading  edge  of  the 
crack  (Figure  13).  The  minor  axis  of  the  disbond  is  along  a  direction  that  is  at  right 
angles  to  the  crack  plane  and  at  the  crack  centerline.  The  ratio  of  minor  to  major  axis  of 
the  disbond  depends  on  a  number  of  parameters  such  as  the  type  of  adhesive  (ductile 
versus  brittle),  adhesive  thickness  and  the  properties  of  bonded  (metal  and  composite) 
materials. 


Location  of  Crack  Tips  in  Cracked  Sheet 


Elliptical  Shape  Disbond  in  Adhesive 

Figure  13.  Schematic  Representation  of  Adhesive  Disbond  Under  Cyclic  Loading 


3.6.  Development  of  Computer  Program  to  Calculate  Normalized  SIF’s 

A  computer  program  was  written  to  compute  shear  stresses  in  the  adhesive  and  stress 
intensity  factors  in  the  cracked  sheet  due  to  applied  loading  and  thermal  mismatch 
between  the  metallic  sheet  and  composite  repair  patch.  The  stress  intensity  factors  can  be 
used  in  making  crack  growth  and  residual  strength  predictions. 

The  computer  program  is  written  to  account  for  the  elliptical  disbond  in  the  adhesive  with 
major  axis  along  the  crack  plane  and  leading  edge  of  the  disbond  coinciding  with  the 
crack  tip.  The  user  has  the  option  to  input  the  minor  to  major  axis  ratio  of  the  disbond  b/a 
(Figure  13).  The  recommendation  for  the  size  of  disbond  to  be  used  in  the  analysis  is 
given  in  the  section  on  design  guidelines  and  in  the  AFGROW  [26]  software  program 
further  developed  as  part  of  the  present  research. 

The  normalized  stress  intensity  factor  is  calculated  based  on  a  unit  applied  remote  stress 
to  an  infinitely  wide  center-cracked  panel.  In  AFGROW,  this  normalized  stress  intensity 
factor  is  equivalent  to  the  ratio  of  the  patched  stress  intensity  factor  to  the  unpatched 
stress  intensity  factor. 

The  computer  program  also  accounts  for  the  effect  of  bending  in  the  structure.  The 
bending  correction  discussed  in  Section  3.3  has  been  incorporated  in  the  program.  The 
computer  program  has  the  following  inputs; 
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Cracked  Sheet  Parameters: 


1)  Elastic  Modulus,  E, 

2)  Poisson  ratio,  v, 

3)  Coefficient  of  linear  expansion,  a, 

4)  Thickness,  hj 

Repair  Patch  Parameters: 

1)  Elastic  modulus  of  the  repair  patch  along  the  crack  plane  direction,  E^ 

2)  Elastic  modulus  of  the  repair  patch  along  loading  direction  at  right  angles  to  crack 
plane,  Ey 

3)  Poisson  ratio, 

4)  Shear  Modulus  G,(y 

5)  Coefficient  of  linear  expansion  of  patch  material  tty  along  loading  direction 

6)  Patch  width  W 

7)  Patch  Length  L 

8)  Patch  Thickness 

Adhesive  Parameters 

1)  Thickness,  ‘h^. 

2)  Shear  Modulus,  ‘pg. 

3)  Ratio  of  minor  to  major  axis  of  debond,  ‘deb’ 

Crack  Parameters 

1)  Half  crack  length,  ‘a’ 

2)  Increment  in  half  crack  length,  ‘da’,  for  computing  stress  intensity  factors 

3)  Number  of  runs,  ‘nab’ 

Temperature  Parameters 

1)  Temperature  difference,  AT,  which  causes  thermal  mismatch 


3.7.  Implementation  of  Composite  Repair/Analysis  into  AFGROW 

This  section  describes  the  implementation  of  a  preliminary  composite  bonded  repair 
design  module  into  the  AFGROW  computer  program.  This  feature  has  been  integrated 
with  the  existing  life  prediction  capability  so  that  the  effect  of  a  repair  patch  may  be 
included  in  crack  growth  analysis.  A  limited  “knowledge  base”  has  also  been  built  into 
this  system  to  aid  the  user  in  selecting  an  initial  patch  material  and  ply  la5Tip.  The 
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following  sections  contain  a  description  of  the  features,  capabilities  and  limitations  of  this 
module. 


3.8.  Repair  Design  Considerations 

A  number  of  factors  need  to  be  considered  in  designing  composite  repair  patches  for 
metallic  structures.  The  following  factors  are  taken  into  account  in  AFGROW; 

1 .  Repair  patch  material  and  ply  thickness. 

2.  Repair  patch  layup. 

3 .  Substructure  material  properties  and  thickness. 

4.  Effect  of  thermal  stresses. 

5.  Effect  of  bending. 

6.  Adhesive  properties  and  thickness. 

7.  Local  disbonding  of  adhesive  around  crack. 

While  there  are  several  other  important  considerations  such  as  repair  location,  repair 
geometry  and  dimensions,  and  edge  taper,  these  are  outside  the  scope  of  the  current 
version  of  AFGROW  (Version  3.85.6).  The  program  is  intended  to  provide  preliminary 
solutions  to  repair  design  problems  and  to  assess  the  effect  of  various  parameters  on  the 
life  of  repaired  substructures.  More  detailed  analyses  would  require  the  use  of  tools 
capable  of  creating  and  analyzing  detailed  models  of  the  repaired  structure.  Such  tools 
may  then  be  used  in  combination  with  AFGROW  for  life  prediction  purposes. 


3.9.  Program  Implementation 

The  bonded  repair  module  was  designed  to  work  with  the  existing  crack  growth  analysis 
capability  so  that  a  repair  can  be  applied  to  the  current  analysis  model  at  any  time.  The 
life  prediction  routine  includes  the  effect  of  the  repair  in  all  subsequent  analyses  until  the 
repair  is  deleted  by  the  user.  However,  repairs  are  not  applicable  to  all  the  crack 
configurations  offered  in  AFGROW.  The  range  of  applicability  is  outlined  in  the  section 
below. 


3.9.1.  Repair  Application  Bounds 

Repairs  are  applicable  to  all  the  through-crack  solutions  in  AFGROW  with  the  exception 
of  the  solid  circular  shaft,  pipe,  and  disc  solutions.  They  may  not  be  applied  to  any  of  the 
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part-through  crack  configurations.  In  addition,  certain  other  limitations  are  placed  on  the 
dimensions  of  the  model  to  be  repaired.  These  are: 

1 .  Minimum  thickness  of  plate  =  0.04” 

2.  Maximum  thickness  of  plate  =  0. 125” 

3 .  Maximum  crack  length  =  2.00” 

The  minimum  thickness  and  maximum  crack  length  limitations  are  due  to  the  range  of 
applicability  of  the  repair  design  solution  module.  The  maximum  thickness  limitation  is 
imposed  because  the  current  release  of  AFGROW  (Version  3.85.6)  can  not  model  the 
variation  in  crack  length  (through  the  thickness)  that  occurs  when  repairs  are  applied  to 
thicker  specimens. 

AFGROW  prohibits  the  application  of  a  repair  if  the  model  configuration  and  dimensions 
are  not  within  the  acceptable  limits. 


3,9,2.  Design  Parameters  and  Controls 

It  is  recommended  to  set  up  all  the  model  parameters  and  properties  before  designing  a 
repair.  AFGROW  uses  these  parameters  and  the  repair  patch  properties  to  compute  a  set 
of  “beta  correction  factors”  which  are  used  to  correct  the  stress  intensity  factors  during 
crack  propagation  to  account  for  the  presence  of  the  repair  patch.  These  factors  are  stored 
in  a  tabular  form  as  a  function  of  crack  length.  If  the  model  geometry,  dimensions  or 
material  properties  are  changed  after  a  repair  has  been  applied,  AFGROW  automatically 
deletes  the  repair  as  the  beta  correction  factors  are  no  longer  valid.  In  such  a  case,  the 
repair  must  be  redesigned. 

Spectrum  information  is  used  by  AFGROW  to  automatically  design  a  repair.  If  the 
“auto-design”  feature  is  to  be  used,  then  the  correct  spectrum  must  be  opened  and  the 
stress  multiplication  factor  must  be  specified  prior  to  using  this  feature.  The  program 
determines  the  maximum  stress  level  from  the  spectrum  and  uses  this  information  to 
choose  the  most  appropriate  repair  material  for  the  application. 

The  repair  design  module  is  invoked  by  clicking  on  the  “Design”  option  in  the  “Repair” 
pull-down  menu.  The  repair  design  dialog,  shovm  in  Figure  14,  appears  on  the  screen. 
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Figure  14.  The  Composite  Patch  Design  Dialog 

This  dialog  contains  various  controls  that  can  be  used  to  change  the  repair  patch 
parameters.  These  parameters  and  the  controls  that  are  used  to  modify  them  are 
explained  in  the  folio-wing  sections. 


3.9.2.I.  Material  Properties 

The  current  patch  material  is  displayed  on  the  upper  left  portion  of  the  dialog  box.  A  new 
material  may  be  specified  by  clicking  the  “Change  Material”  button.  This  causes  a  file 
browser  containing  patch  material  data  files  to  pop  up.  This  release  of  AFGROW 
includes  one  such  file,  called  patch .  pch,  which  contains  material  properties  for  boron, 
graphite,  and  GLARE.  The  user  may  then  select  one  of  these  materials  to  be  the  new 
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patch  material.  If  a  material  other  than  these  is  to  be  used,  data  for  that  material  must  be 
present  in  patch .  pch  or  any  other  * .  pch  file. 

The  material  properties  are  referred  to  on  the  dialog  by  means  of  abbreviations,  as  seen  in 
Figure  14.  The  list  of  properties  is: 

EL :  Y oung  ’  s  modulus  in  longitudinal  direction 

ET :  Y oung  ’  s  modulus  in  transverse  direction 

GLT:  Shear  modulus 

PNU:  Poisson’s  ratio 

ALFL;  Coefficient  of  thermal  expansion  in  longitudinal  direction 

ALFT :  Coefficient  of  thermal  expansion  in  transverse  direction 


It  must  be  noted  that  the  program  does  not  allow  the  user  to  alter  the  values  of  the 
material  properties  directly.  This  is  because  the  repair  module  relies  on  the  material 
being  orthotropic;  unexpected  results  may  occur  if  the  properties  do  not  satisfy  the 
conditions  of  orthotropy.  Since  most  repair  applications  utilize  boron,  graphite  or 
GLARE,  the  material  datafile  included  with  AFGROW  is  sufficient  in  most  cases.  Other 
materials  must  be  used  only  after  ensuring  that  they  are  orthotropic. 

A  default  ply  thickness  has  been  set  for  each  of  these  materials.  This  value  is  editable, 
and  may  be  set  to  any  value.  The  other  editable  value  is  AT,  which  is  used  to  account  for 
thermal  stresses  due  to  differences  in  the  curing  temperature  of  the  patch  and  the 
temperature  of  the  plate. 


3.9.2.2.  Patch  Dimensions 

The  patch  dimensions  (width  and  length)  are  used  only  for  the  purposes  of  calculating  the 
effect  of  bending.  The  beta  correction  factors  are  based  on  the  assumption  of  an  infinite 
patch;  the  patch  dimensions  do  not  affect  them  imless  bending  is  included  in  the  solution. 


3.9.2.3.  Ply  Layup 

This  release  of  AFGROW  can  model  patches  with  up  to  32  plies.  The  user  may  specify 
the  orientation  of  each  ply  by  selecting  the  ply  from  the  scrollbox  and  entering  its 
orientation  in  the  textbox  above  it.  If  a  symmetric  patch  with  an  even  number  of  plies  is 
to  be  modeled,  then  it  is  sufficient  to  enter  the  orientations  for  half  the  plies;  AFGROW 
will  automatically  mirror  the  layup  to  construct  the  other  half 
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It  is  recommended  to  have  70%  0°  plies,  20%  ±45°  plies,  and  10%  90  plies.  However, 
in  most  cases  it  is  not  possible  to  achieve  this  distribution  and  the  closest  distribution 
must  be  used.  Also,  it  is  desirable  to  place  the  ±45°  plies  on  the  surface  as  this  will 
reduce  the  tendency  of  the  patch  to  delaminate  at  the  edges. 


3.9.2.4.  Patch  Type 

AFGROW  allows  the  user  to  specify  certain  other  options  that  are  important  from  a 
design  consideration.  These  are  controlled  by  means  of  toggle  buttons  and  include: 

Symmetric  patch 

Specifying  this  option  automatically  mirrors  the  specified  layup  to  create  a  symmetric 
patch.  Consideration  must  be  given  to  maintaining  balanced  laminate  symmetry  to  avoid 
warpage  in  the  precure  of  the  patch. 

Double-sided  patch 

This  option  applies  the  patch  to  both  sides  of  the  metallic  substructure.  Selection  of  this 
option  may  depend  on  factors  such  as  location  and  accessibility  of  the  repair. 

An  assumption  is  made  in  AFGROW  that  the  effect  of  bending  is  negligible  when 
patches  are  applied  to  both  sides  of  the  plate.  The  third  toggle  button,  which  controls 
whether  bending  is  considered  or  not,  is  then  automatically  selected  and  “greyed  out,”  or 
made  inaccessible,  to  the  user. 

No  Bending 

The  user  may  select  this  button  to  either  include  or  neglect  the  effect  of  bending, 
depending  on  whether  bending  due  to  the  presence  of  the  patch  is  deemed  to  be 
significant  at  the  repair  location  or  not. 

It  may  be  turned  on  to  assess  the  contribution  of  bending  to  the  total  stresses  at  the  repair 
location,  and  to  study  the  effect  of  bending  on  the  fatigue  life  of  the  patched  structure. 

It  must  be  noted  that  “bending”  in  this  case  refers  only  to  the  bending  induced  by  the 
repair  patch  on  the  plate. 
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3.9.2.5.  Adhesive  Properties 


The  user  may  specify  the  adhesive  that  is  used  in  the  repair  by  entering  the  name,  shear 
modulus  (GXY),  and  thickness  of  the  adhesive.  In  addition,  a  variable  called  DEB, 
which  characterizes  the  extent  of  the  disbond  in  the  adhesive  in  the  region  of  the  crack, 
must  be  specified.  The  disbond  is  assumed  to  be  elliptical  in  shape,  with  the  crack  lying 
along  the  major  axis  of  the  ellipse.  DEB  is  the  aspect  ratio  of  this  ellipse.  The  value  of 
DEB  must  be  between  0  and  1 . 


3.9.2.6.  The  Patch  Stiffness  Indicator 

Near  the  bottom  of  the  dialog  box,  and  spanning  its  width,  lies  the  patch  stiffiiess 
indicator.  A  diamond-shaped  marker  indicates  the  equivalent  stiffness  of  the  current 
patch  as  a  percentage  of  the  stiffness  of  the  plate  to  which  the  repair  is  to  be  applied. 

It  is  desirable  for  a  composite  repair  patch  to  have  an  equivalent  stiffiiess  between  100% 
and  120%  of  the  stiffness  of  the  metallic  substructure.  In  AFGROW,  this  region  is 
considered  “ideal,”  with  the  70%  -  100%  and  120%  -  150%  ranges  being  allowable,  and 
any  value  outside  these  ranges  considered  unacceptable.  If  the  stiffness  indicator  shows  a 
stiffness  that  is  not  acceptable,  the  program  will  not  allow  that  design  to  be  applied. 

Every  time  any  of  the  parameters  are  changed,  the  equivalent  properties  of  the  patch  are 
recomputed  based  on  laminate  theory,  and  the  stiffness  indicator  moves  to  the  position 
corresponding  to  the  new  equivalent  stiffiiess  of  the  patch. 

The  patch  stiffiiess  indicator  is  intended  as  an  aid  to  the  user  in  monitoring  the  relative 
stiffnesses  of  the  patch  and  the  plate. 


3.9.2.7.  The  Auto  Design  Feature 

Using  this  feature  selects  the  most  suitable  patch  configmation  for  the  current  problem 
definition  and  updates  all  the  controls  on  the  repair  dialog  box  to  reflect  the  new  values. 
Auto  Design  selects  the  patch  material  based  on  the  maximum  stress  in  the  spectrum 
(boron  for  25  KSI  and  above,  graphite  for  20-25  KSI  and  GLARE  for  below  20  KSI)  and 
the  thickness  of  the  metal  plate  (GLARE  for  t  <  1/16”,  graphite  or  boron  for  t  >  1/16”). 
Auto  Design  then  determines  an  optimum  number  of  plies  that  satisfies  the  stiffness 
requirement  while  maintaining  laminate  symmetry  and  a  reasonable  level  of  orthotropy 
(%  of  ±45°  and  90°  plies  as  outlined  in  section  3.10.2.5). 
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3.9.2.8.  Other  Controls 


The  “Critical  SIF”  dialog  allows  the  user  to  control  whether  the  critical  stress  intensity 
factor  will  be  based  on  the  patched  structure  or  the  unpatched  structure,  by  selecting  the 
appropriate  toggle  button.  While  this  has  no  bearing  on  the  patch  design  itself,  it  is  used 
during  fatigue  crack  propagation  to  determine  the  failure  criteria. 

“Show”  buttons  are  used  to  select  the  plots  that  are  to  be  displayed  after  the  beta 
correction  factors  are  computed.  These  include  “Beta  Correction  Factor  plots  and  Stress 
Contour  plots  of  the  stresses  in  the  patch  and  the  adhesive. 


3.9.3.  Patch  Design 


After  all  the  design  parameters  have  been  specified,  the  “Apply”  button  may  be  clicked  to 
compute  the  beta  correction  factors.  At  this  point  the  main  AFGROW  window  displays  a 
color-coded  schematic  of  the  ply  orientations,  and  a  list  of  the  equivalent  properties  of  the 
composite  patch,  which  are  computed  using  laminate  theory.  Figure  15. 


Matenal  Name:  BORONCPOXY 

EX-  1.65671  e+07 

'EY  -  9,S4469e+06 

QXY  -  3.308078+06 

THICK  »  0.0572  ■ 

NUXV'  -  0.314233 
MJ-TX-  0.181037  , 

ALFL  -  2.76384e-0B 
ALFT-  3.861 02e-D6 


Beta  Correction  Factors 
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Figure  15.  Ply  Orientations  and  Equivalent  Patch  Properties 


The  beta  correction  factors  are  plotted  on  this  screen  as  soon  as  they  are  computed.  A 
bigger  window  that  allows  the  user  to  overlay  several  beta  correction  factor  plots  is  also 
displayed  if  the  option  is  selected.  Figure  16. 
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Figure  16.  Beta  Correction  Factor  Plot  Window 


The  program  also  displays  contours  of  the  shear  stresses  in  the  adhesive  and  the  stresses 
in  the  patch  in  a  separate  window,  Figure  17: 
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Figure  17.  Stress  Contour  Window 


This  window  includes  toggle  buttons  to  display  either  the  adhesive  or  patch  stresses. 
Currently,  this  release  of  AFGROW  does  not  include  the  capability  to  display  skin 
stresses  and  the  toggle  button  is  made  inactive. 

The  slider  bar  at  the  lower  right-hand  comer  controls  the  crack  length.  As  the  stresses  in 
the  patch  and  the  adhesive  redistribute  continuously  as  the  crack  grows,  this  effect  may 
be  simulated  by  dragging  the  slider  bar  to  the  desired  crack  length. 

It  must  be  noted  that  the  size  of  the  plot  area  and  the  representation  of  the  crack  are  kept 
constant  when  the  crack  length  is  changed  using  the  slider  bar.  A  display  of  the  current 
coordinates  and  the  stress  at  that  point  are  provided  to  the  left  of  the  plot  when  the  cursor 
is  brought  into  the  plot  window.  In  addition,  the  stresses  shown  are  for  the  center- 
cracked  case  and  are  not  intended  to  reflect  the  actual  stresses  present  for  other 
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geometries.  However,  it  is  assumed  that  the  Beta  correction  factors  will  provide  a 
reasonable  approximation  of  stress  intensity  for  other  geometries,  especially  if  the  cracks 
have  grown  a  reasonable  distance  from  the  geometric  configuration  (notch,  hole,  etc.) 

When  a  satisfactory  patch  design  is  achieved,  the  patch  design  dialog  may  be  dismissed. 
At  this  point,  a  representation  of  the  patch  on  the  plate  is  displayed  in  the  main  window 
as  shown  in  Figure  1 8. 


<=  AFGftOWv3.83Beta:Untmed  »  Q 


Figure  18.  Main  Window  Showing  Applied  Patch 

All  subsequent  crack  growth  analyses  will  include  the  effect  of  the  patch  imtil  it  is  deleted 
using  the  “Dismiss”  option  in  the  Repair  menu. 
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The  patch  is  automatically  deleted  when  the  model  configuration,  dimensions,  or  material 
properties  are  changed.  This  is  because  the  stored  beta  factors  are  design  specific  and  no 
longer  valid. 

The  effect  of  different  patch  designs  on  fatigue  crack  growth  life  may  be  observed  by 
overlaying  the  crack  growth  plots.  Other  features  of  the  code,  such  as  retardation  models, 
residual  stress  capability,  and  extreme  environment  effects  may  be  used  simultaneously 
with  the  composite  repair  capability. 
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4. 


WEIGHT  FUNCTION  DEVELOPMENT 


Durability  and  strength  evaluation  of  notched  and  cracked  structural  elements  requires  the 
calculation  of  stress  intensity  factors  for  cracks  located  in  regions  characterized  by 
complex  stress  fields.  This  is  particularly  true  for  cracks  emanating  from  notches  and 
other  stress  concentration  regions  that  might  be  frequently  foimd  in  mechanical  and 
structural  components.  In  the  case  of  engine  components,  complex  stress  distributions 
are  often  due  to  temperature  effects.  The  existing  stress  intensity  factor  handbook 
solutions  are  not  sufficient  in  such  cases  due  to  the  fact  that  most  of  them  have  been 
derived  for  simple  geometry  and  loading  configurations.  The  variety  of  notch  and  crack 
configurations  and  the  complexity  of  stress  fields  occurring  in  engineering  practice 
require  more  efficient  tools  for  calculating  stress  intensity  factors  than  the  large  but 
nevertheless  limited  number  of  currently  available  ready  made  solutions  obtained  for  a 
few  specific  geometry  and  load  combinations. 

Therefore,  methods  were  developed  for  the  efficient  calculation  of  stress  intensity  factors 
for  cracks  in  complex  stress  fields  by  using  the  weight  function  approach. 


4.1.  Technical  Background 

Most  of  the  existing  methods  of  calculating  stress  intensity  factors  require  a  separate 
analysis  for  each  load  and  geometry  configuration.  The  weight  fimction  method 
developed  by  Bueckner  [27]  and  Rice  [28]  simplifies  considerably  the  determination  of 
stress  intensity  factors.  If  the  weight  function  is  known  for  a  given  cracked  body,  the 
stress  intensity  factor  due  to  any  load  system  applied  to  the  body  can  be  determined  by 
using  the  same  weight  function.  Therefore,  there  is  no  need  to  derive  a  ready  made  stress 
intensity  factor  for  each  load  system  and  resulting  internal  stress  distribution.  The  stress 
intensity  factor  can  be  obtained  by  multiplying  the  weight  function,  m(x,a),  and  the 
internal  stress  distribution  a(x)  in  the  prospective  crack  plane  and  integrating  the  product 
along  the  crack  length  ‘a’. 

The  success  of  the  weight  function  technique  for  calculating  stress  intensity  factors  lies  in 
the  possibility  of  using  superposition.  It  can  be  shown  [29]  that  the  stress  intensity  factor 
for  a  crack  body  (Figure  19)  subjected  to  the  external  loading,  S,  is  the  same  as  the  stress 
intensity  factor  in  a  geometrically  identical  body  with  the  local  stress  field,  ci(x),  applied 
to  the  crack  faces.  The  local  stress  field,  a(x),  in  the  prospective  crack  plane  is  due  to  the 
external  load,  S,  and  it  is  determined  for  the  uncracked  body  by  ignoring  the  presence  of 
the  crack. 
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(a)  Local  stress  field  in  an  uncracked  body 
subjected  to  external  loading,  S 


(b)  Geometrically  identical  body  with  local 
tractions,  a(x)  applied  to  crack  faces 


Figure  19.  Nomenclature  and  Concept  of  Superposition 
(used  with  permission  from  G.  Glinka) 


The  vinique  feature  of  the  weight  function  method  is  that  once  the  weight  function  for  a 
particular  cracked  body  is  determined,  the  stress  intensity  factor  for  any  loading  system 
applied  to  the  body  can  be  calculated  by  simple  integration  of  the  product  of  the  stress 
field,  c(x),  and  the  weight  function,  m(x,a). 


K  =  ni(x,  a)  dx  (34) 

The  idea  of  using  weight  functions  is  illustrated  in  Figures  20a  and  20b,  where  the 
through  thickness  crack  is  used  as  an  example.  The  weight  function,  m(x,a),  can  be 
interpreted  as  the  stress  intensity  factor  that  results  from  a  pair  of  splitting  forces  applied 
to  the  crack  face  at  position  x.  Since  the  stress  intensity  factors  are  linearly  dependent  on 
the  applied  loads,  the  contributions  from  multiple  splitting  forces  applied  along  the  crack 
surface  can  be  superposed  and  the  resultant  stress  intensity  factor  can  be  calculated  as  the 
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sum  of  all  individual  load  contributions.  This  results  in  the  integral  (34)  of  the  product  of 
the  weight  function,  m(x,a),  and  the  stress  field,  a(x),  for  the  case  of  continuously 
distributed  stress  field. 

The  Weight  Function  Approach  to 
Calculation  of  Stress  Intensity  Factors 

Weight  Function  -  stress  intensity  factor,  K,  for  a 
pair  of  splitting  forces  applied  to 
crack  surfaces 


For  a  through  crack 
in  infinite  plate: 

/  1  la  +  x 

m^(x,a)  =  -—  - 

WTta  \a-x 

■Jna  Va  +  x 


Stress  intensity  factors  for  crack  crack  tips  A  and 
B  respectively: 

f'  /  m  P  1(2  +  X 

VJca  \a-x 

Kg  =  mg  {x,  a,  P)  =  pHZ 
Vita  Va  +  x: 


Figure  20a.  Weight  Function  and  Stress  Intensity  Factors  for  a  Through-Crack  in  an  Infinite  Plate 

(used  with  permission  from  G.Glinka) 
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The  Weight  Function  Approach  to 
Calculation  of  Stress  Intensity  Factors 


Tn  liie  case  of  two  forces^  by  using  superposition 


■yjiui  ^  fli  -  jCj  vTia  ^  a  - 

In  the  case  of  multiple  forces 

Ka  =  rn{a,  jtj,  +  w(fl,  X2 ,  .  A-m{a,  Pfi 

For  the  stress  distribution  a(x) 

a 

=  \^ix)m^ix,a)dx 


Kg  -  ja{x)mif{x^a)dx 
-a 


Figure  20b.  Weight  Function  and  Stress  Intensity  Factors  for  a  Through-Crack  in  an  Infinite  Plate 
Subjected  to  multiple  Loads  and/or  Continuously  Distributed  Stress. 

(used  with  permission  from  G.  Glinka) 
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4.2.  Universal  Weight  Functions  for  One-Dimensional  Cracks 

The  weight  function  is  dependent  on  the  geometry  only  and  in  principle  should  be 
derived  individually  for  each  geometrical  configuration.  However,  Glinka  and  Shen  [30] 
have  found  that  one  general  weight  function  expression  can  be  used  to  approximate 
weight  functions  for  a  variety  of  geometrical  configurations  of  cracked  bodies  with  one 
dimensional  cracks  of  Mode  I: 


The  system  of  coordinates  and  the  notation  for  through  internal  and  edge  cracks  are  given 


in  Figure  21. 
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Figure  21.  Single  Edge  Crack  and  Central  Through-Crack  Under  Symmetrical  Loading 

In  order  to  determine  the  weight  function,  m(x,a),  of  eq.  (35)  for  a  particular  cracked 
body  it  is  sufficient  to  determine  [31]  the  three  parameters  M„  M2,  and  Mj.  Because  the 
weight  function  of  eq.  (35)  is  the  same  for  all  cracks  then  the  same  integration  routine  can 
be  used  for  calculating  stress  intensity  factors  from  eq.  (34).  Moreover,  only  a  limited 
number  of  generic  weight  functions  is  needed  to  enable  the  determination  of  stress 
intensity  factors  for  a  large  number  of  load  and  geometry  configurations. 
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4.3.  Universal  Weight  Functions  for  Two-Dimensional  Part-Through  Surface  and 
Corner  Cracks 


In  the  case  of  2-D  cracks  such  as  semielliptical  and  comer  surface  cracks  in  plates  and 
cylinders,  the  stress  intensity  factor  changes  along  the  crack  front.  However,  in  most 
practical  cases  the  deepest  point  A  (Figure  22)  and  the  surface  point  B  are  associated  with 
the  highest  and  the  lowest  value  of  the  stress  intensity  factor  along  the  crack  front. 
Therefore,  the  weight  functions  for  the  points  A  and  B  have  been  derived  [32] 
analogously  to  the  rmiversal  weight  function  (35). 


For  point  A  (Figure  22) 
2 


mA(x,a,a/c,a/t)  =  ' 


^J2n{a.-x) 
•  For  point  B  (Figure  22) 


1  +  M 


lA 


2 

1 

f.  x\ 

1-- 

^2A 

1-- 

+  M3A 

1  — 

V  aj 

l  aj 

L  aj 

mB(x,a,a/c,a/t)  = 


1  +  M 


IB 


fx^ 

2  1 

fx^ 

I 

fx] 

^2b| 

^3B 

UJ 

(36) 


(37) 


The  weight  functions  mpi[x,a)  and  mB(x,a)  given  above  corresponding  to  the  deepest  and 
the  surface  points  A  and  B,  respectively,  have  been  derived  for  one-dimensional  stress 
fields  (Figure  22),  dependent  on  one  variable ,  x,  only. 
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Figure  22.  Semielliptica!  Surface  Crack  in  an  Infinite  Plate  of  Finite  Thickness,  t 


4.4.  Sequence  of  Steps  for  Calculating  Stress  Intensity  Factors  Using  Weight 

Functions 

In  order  to  calculate  stress  intensity  factors  using  the  weight  function  technique  the 
following  tasks  need  to  be  carried  out: 


•  Determine  stress  distribution  a(x)  in  the  prospective  crack  plane  using  linear  elastic 
analysis  of  uncracked  body  (Figure  19a),  i.e.,  perform  the  stress  analysis  ignoring  the 
crack  and  determine  the  stress  distribution  a(x)  =  Gq  f(S,x);  where:  Go  -  nominal 
stress,  X  -  coordinate,  S  -  external  load 


•  Apply  the  “uncracked”  stress  distribution,  cr(x),  to  the  crack  surfaces  (Figure  19b)  as 
tractions 
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•  Choose  appropriate  generic  weight  function  (i.e.,  choose  appropriate  M,,  M2,  and  M3 
parameters). 

•  Integrate  the  product  of  the  stress  function  a(x)  and  the  weight  function  m(x,a)  over 
the  entire  crack  length  or  crack  surface,  eq.  (34). 


4.5.  Determination  of  Weight  Functions 

The  derivation  of  the  weight  function  of  eq.  (35)  for  a  one-dimensional  crack  is 
essentially  reduced  to  the  determination  of  parameters  M„  Mj,  and  M3.  There  are  several 
methods  available  depending  on  the  amount  of  information  available  as  described  in  [31]. 
The  most  frequently  used  method  based  on  two  reference  stress  intensity  factors  is 
described  below. 

Supposing  that  two  reference  stress  intensity  solutions  K^,  and  Kj2  for  two  different  stress 
fields  cr^i(x)  and  cr,2(x),  respectively,  are  known,  a  set  of  two  equations  can  be  written  by 
using  eqs.  (34)  and  (35).  The  third  equation  necessary  to  determine  parameters  M„  Mj, 
and  M3  can  be  formulated  based  on  the  knowledge  of  the  crack  surface  slope  [31,32]  at 
the  crack  center  (central  through  cracks)  or  from  the  crack  surface  curvature  at  the  crack 
mouth  (edge  cracks).  The  three  pieces  of  information  result  in  three  independent 
equations  which  can  be  used  for  the  determination  of  unknown  parameters  M„  Mj,  and 

M3; 

•  Central  Through  Cracks  Under  Symmetric  Loading  (Figure  21a) 


Krl  =  J[<^rl(x)- 


2  L . f.  xv  ,,  r,  x^  r,  xv 


Ima  -  X 


1  +  M,  1 —  +M2  1--  +M3  1--  dx  (38) 


Kr2  =  J['^r2(x)- 


2  f  xV  f  f  xM 

.  ,  1  +  M,  1 - +M2  1--  +M3  1--  dx  (39) 

Z7t(a-x)  I  aj  n  aj  T  a 


50 


Edge  Cracks 


K„  =  fcr.,(x)- 


0  {  v\2  (  f  X12 

■  ^  1  +  Mj  1  —  +  M2  1  —  +  M3  1  —  dx 

Z7i(a~x)  V  V  V 


K.  =  fa.(x)^=^  1  +  M.  1-iJ+Mil-i  +M3  1-fJ  dx  (42) 


iTUla-x 


a'  2 

5x^  1  J2n(a.-x] 


l  +  M,(l--r +Mjl--1  +M,(l--1^ 


The  two  reference  cases  can  be  obtained  from  the  literature  or  they  can  be  derived  using 
finite  element  method.  Most  often  available  are  the  solutions  for  the  uniform  and  linear 
stress  distributions  induced  by  simple  tension  or  bending  load  and  those  cases  were 
predominantly  used  to  derive  the  weight  fimctions  listed  below. 

•  Semielliptical  Surface  and  Corner  Cracks 

The  set  of  equations  necessary  for  deriving  parameters  for  the  weight  function 
corresponding  to  the  deepest  point  A  on  the  crack  front  (Figure  21)  was  found  to  be 
identical  to  the  set  of  equations  derived  for  the  edge  crack. 


Kf,  =  fa„(x)- 


y  yA2  f  f  XU 

I  +  Mja  1--  +M2A  1--  +M3A  1--  dx  (44) 

V  a)  V  ay  V  ay 


^ 


In  the  case  of  the  surface  point  B  the  first  two  equations  were  the  same  as  previously  and 
the  additional  equation  was  derived  [32]  by  satisfying  the  condition  that  the  weight 
function  must  vanish  for  x=0. 
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dx  (47) 
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dx 


(48) 


0  - 1  +  Mjb  +  M2B  +  Mjb  (49) 

The  three  unknown  parameters  Mj  can  be  determined  by  simultaneously  solving  one  of 
the  set  of  equations  presented  above. 


4.6.  Numerical  Integration  of  the  Weight  Function  and  Calculation  of  Stress 
Intensity  Factors 

The  calculation  of  stress  intensity  factors  from  the  weight  function  requires  integration  of 
of  the  product  “s(x)m(x,a)”  along  the  crack  length  according  to  eq.  (34).  The  weight 
function  can  always  be  written  in  the  general  form  of  eq.(36)  or  eq.(37).  However,  the 
stress  distribution  “s(x)”  can  take  any  form  depending  on  the  problem  of  interest.  If  the 
stress  distribution  is  given  in  the  form  of  a  mathematical  expression,  analytical 
integration  can  be  performed  and  closed  form  integrals  of  eq.  (1)  are  sometimes  feasible. 
However,  very  often  the  stress  distribution  “s(x)”  is  obtained  from  finite  element 
calculations  and  the  results  are  given  as  a  series  of  stress  values  corresponding  to  a  range 
of  point  of  coordinate  “x.”  Therefore,  a  numerical  integration  technique  is  needed  for  the 
integration  of  equation  (34)  and  the  calculation  of  stress  intensity  factors.  Two  methods 
of  efficient  integration  of  eq.  (34)  are  described  below. 


4.6.1.  Integration  Using  Centroids  of  Areas  Under  the  Weight  Function  Curve 

The  integration  method  using  the  area  centroids  is  based  on  the  following  theorem: 

If  m(x,a)  and  C7(x)  are  monotonic  and  linear  function,  respectively,  and  both  depend  on 
variable  x  only,  (Figure  23),  then  the  integral  (34)  can  be  calculated  from  expression  (47), 
representing  the  product  of  the  area  S  under  the  curve  m(x,a)  and  the  value  of  the  function 
a(X)  corresponding  to  the  coordinate  x=X  of  the  centroid  C. 

K  =  S*a(x)  (50) 

The  weight  functions,  m(x,a),  are  monotonic  and  nonlinear.  The  stress  functions,  a(x). 
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are  usually  nonlinear  as  well.  Therefore,  in  order  to  apply  the  theorem  above  to  the 
integral  (34),  the  integration  interval  is  divided  into  “n”  subintervals  in  such  a  way  that, 
the  stress  function  cy(x)  is  approximated  by  the  secant  line  drawn  between  the  end  points 
of  each  subinterval  Figure  22.  Thus,  the  stress  function  a(x)  over  the  subinterval  “I,” 
may  be  written  in  the  form  of  eq.  (51), 

ai(x)  =  AjX  +  Bj  for  (xj -l)<x<Xi  (51) 

where : 

^  ^  g(xi)-g(xi  j  =a(xi)-AiXi  (51a) 

'  Xi-Xi_, 

After  substitution  of  expression  (51)  into  eq.  (34)  and  summation  over  all  subintervals  the 
following  expression  for  the  stress  intensity  factor  can  be  derived. 

K  =  ^  |(AiX  +  Bi)m(x,a)dx  (52) 

i=I  x,., 

Each  integral  in  eq.  (52)  can  be  computed  by  using  the  simplified  (Figure  23)  integration 
method  given  in  the  form  of  eq.  (50).  Thus,  the  stress  intensity  factor,  K,  can  be  finally 
written  in  the  form  of  expression  (53). 


K  =  ^Si*a(Xi),  where  i  =  1,2,.. .n  (53) 

i 

In  order  to  calculate  the  stress  intensity  factor  given  in  the  form  of  eq.  (53),  it  is  necessary 
to  calculate  the  areas  Sj  under  the  weight  function  curve  m(x,a),  and  the  coordinates  of 
their  centroids,  Xj.  Both  the  areas  S;  and  the  centroid  coordinates  Xj  for  each  sub-interval 
can  be  calculated  once  in  a  general  form  based  on  the  generalized  weight  function  (35). 
However,  they  appear  to  be  too  lengthy  for  an  easy  hand  calculation.  F ortunately ,  further 
simplification  of  the  integration  routine  is  possible  due  to  the  fact  that  the  weight 
functions  are  smooth  within  their  ranges  of  integration.  Therefore  the  procedure  can  be 
reversed  (Figure  25)  by  calculating  first  the  areas  Sj*  under  the  stress  function  a(x)  and 
the  coordinates  X'  of  their  centroid.  Then  the  appropriate  values  of  the  weight  function 
m(X*,a)  can  be  calculated  from  expression  (35).  It  is  worth  noting  that  in  the  case  of  the 
piece-wise  approximation  of  the  stress  function,  a(x),  the  areas  S;*  and  the  coordinates  of 
their  centroids,  X;’,  can  be  easily  calculated  from  relations  (54)  and  (55),  respectively. 
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(54) 


S’  =|Wxi)+a(xj_,)Xxi  -Xi_,) 


X  =x. 


(Xj  -Xi_iX2a(xi_,)+a(xi)] 
3[a(xi)+cj(xi.,)] 


(55) 


Finally  the  stress  intensity  factor,  K,  is  calculated  from  expression  (56). 


K  =  y;s:m(x:,a)  (56) 

i=l 

Thus,  the  numerical  procedure  for  calculating  the  stress  intensity  factor  using  the 
integration  method  described  above  requires  the  calculation  of  appropriate  parameters 
using  equations  (35)  and  (54-56).  The  method  described  above  is  recommended  for  quick 
approximate  calculations  with  the  help  of  a  hand  calculator. 
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Figure  23.  Graphical  Representation  of  Simplified  Integration  of  the  Product  of  Two  One- 
Dimensional  Functions,  a)  Monotonic  Nonlinear  Weight  Function  m(x,a),  b)  Linear  Stress  Function, 

S(x) 
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Figure  25.  Integration  Method  Based  on  the  Subareas  and  Centroids  of  the  Linearized  Piecewise 
Stress  Function,  a)  Weight  Function  ni(x,a),  Linearized  Piecewise  Stress  Function 


4.6.2.  Analytical  Integration  of  the  Linearized  Piecewise  Stress  Distribution  and 
the  Weight  Function 

The  integration  technique  described  in  Section  4.6.1  is  convenient  for  hand  calculator 
calculation  when  the  stress  distribution  can  be  approximated  by  a  few  linear  segments. 
However,  the  segments  adjacent  to  the  crack  tip  can  not  be  large  because  the  weight 
function  tends  to  infinity  near  the  erack  tip  and  if  the  stress  is  the  highest  near  the  crack 
tip  the  method  described  above  might  be  sometimes  inaccurate.  Moreover,  in  the  case  of 
stress  distributions  eharacterized  by  high  gradients  aceurate  approximation  of  the  stress 
distribution  requires  relatively  large  number  of  linear  pieces  and  the  integration  has  to  be 
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carried  out  with  the  help  of  a  computer.  However,  the  computer  integration  routine  can 
also  be  significantly  simplified  because  closed  form  solution  to  the  integral  (52)  can  be 
derived  analytically  for  each  linear  piece  (Figure  25)  of  the  stress  function  s(x). 

The  stress  fimction  s(x)  over  the  linear  segment  “i”  can  be  given  in  the  form  of  the  linear 
equation  (51).  Thus  the  contribution  to  the  stress  intensity  factor  associated  with  the 
stress  segment  “i”  can  be  calculated  from  eq.  (34)  after  substituting  appropriate 
expressions  for  the  stress  and  the  weight  fimction.  The  solutions  given  below  have  been 
derived  for  the  deepest  and  for  the  surfaee  point  of  semielliptical  crack  using  the  weight 
function  (36)  and  (37),  respectively.  The  solution  to  one-dimensional  cracks  is  the  same 
as  for  the  deepest  point  of  semielliptical  cracks. 

•  Deepest  Point  A  (Figure  22) 


K*  =  £(A,x  +  Bi)- 


27ca 


^  x^ 
1-- 


1  +  M 
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1 

,  3“ 
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(.  x) 

(  xV 

1  — 

+  ^2A 
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+  M3A  1  — 

aj 

1  a) 

'  V  a; 

dx 


V  V  a; 

•  Surface  Point  B  (Figure  22) 


K'  =  p  (A,x  +  B,)-^ 

V27TX 


1  +  M 


IB 


^X\2 


+  M,J-|  +M3B 


^x\2 


dx  (58) 


The  closed  form  expressions  resulting  from  the  integration  of  eq.  (57)  and  eq.  (58)  are 
given  below. 


(57) 
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Deepest  Point  A  (Figure  22) 


Kr=j— 


2  faiCCi, +MjACi2  h-MjaCo +M3ACj4) 


+  |3,(C„  +  M,,C„  +  M„C,5  +  M„C„) 


where ; 

ttj  =  Bj  +  aAj  and  P;  =  -aA; 
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Surface  Point  B  (Figure  22) 
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(60) 


where : 

tt;  =  Bj  +  aAj  and  P;  =  -aAj 
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Equations  (59)  and  (60)  can  be  used  for  calculating  stress  intensity  contributions  due  to 
each  linear  piece  of  the  stress  distribution  function  by  substituting  appropriate  values  for 
a,  X;.,,  Xj,  A;  and  Bj.  The  stress  intensity  factor  K  can  be  finally  calculated  as  the  sum  of  all 
contributions  Kj  associated  with  the  linear  pieces  within  the  range  of  0  =  x  =  a. 


(61) 


Thus,  the  integration  can  be  reduced  to  the  substitution  of  appropriate  parameters  into 
equations  (59-60)  and  summation  according  to  eq.  (61).  This  makes  it  possible  to 
develop  very  efficient  numerical  integration  routines,  which  is  important  in  the  case  of 
lengthy  fatigue  crack  growth  analyses. 


4.7.  Implementation  of  Weight  Function  Methods  into  AFGROW 

This  section  describes  the  adaptation  of  the  Weight  Function  Solutions  into  the 
AFGROW  program.  It  also  describes  the  procedure  required  to  apply  a  stress  distribution 
to  a  weight  function  model  configuration  in  the  AFGROW  program. 
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4.7.1.  Selecting  a  Weight  Function  Stress  Intensity  Solution  Model 

To  select  a  Weight  Fvinction  Stress  Intensity  Solution  model;  go  to  the  Insert  menu; 
select  Model;  and  then  select  Weight  Function  Solutions  menu  item.  This  will  cause  a 
Model  Configuration  Dialog  Box  (Figure  26)  to  appear  on  the  screen. 


Figure  26.  Weight  Function  Model  Configuration  Dialog 


This  dialog  consists  of  a  scrolled  window  that  displays  iconic  representations  of  the  nine 
Weight  Function  model  configurations  currently  available  in  the  AFGROW  program.  To 
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the  right  of  each  of  these  icons  is  a  brief  description  of  the  corresponding  model 
configuration.  You  may  select  any  of  the  weight  function  solutions  in  this  dialog  by 
selecting  the  toggle  button  to  the  left  of  the  model  of  interest.  Depressing  any  toggle 
button  causes  any  toggle  button  that  is  already  depressed  to  pop  out,  ensuring  that  only 
one  configuration  is  selected. 

After  selecting  the  desired  configuration,  the  user  must  click  the  DONE  button  at  the 
bottom  of  the  dialog  to  choose  that  configuration  for  further  analysis. 

When  a  Weight  Function  Stress  Intensity  Solution  model  has  been  chosen,  the  title  of  the 
specimen  configuration  in  the  Main  Display  Area  will  include  a  “Weight  Function 
Solution”  statement  (Figure  27). 


Quarter-elliptical  Corner  Crack 
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Figure  27.  AFGROW  Model  Configuration  Display 


4.7.2.  Entering  a  Stress  Distribution 

Each  Weight  Function  Stress  Intensity  Solution  Model  requires  an  “uncracked”  stress 
distribution  along  crack  plane. 

The  AFGROW  program  uses  a  normalized  stress  distribution  in  the  prospective  ‘a’  crack 
plane  (Sy)  to  calculate  Stress  Intensity  Factors  (SIF)  for  the  part-through  cracks  (2-D 
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solutions);  and  a  normalized  stress  distribution  in  prospective  ‘c’  crack  plane  (SJ  to 
calculate  SIF  for  the  through-cracks  (1-D  solutions). 

Note  that  both  solutions  currently  allow  the  input  stress  field  to  vary  in  the  one  direction 
only.  When  a  part-through  crack  transitions  into  a  through-the-thickness  crack,  the  model 
is  automatically  converted  into  the  appropriate  1-D  case,  and  the  applicable  stress 
distribution  is  used  to  continue  the  life  prediction. 

It  is  recommended  that  the  weight  fimction  stress  distribution(s)  should  be  normalized 
according  to  the  stress  value  at  the  crack  origin.  However,  it  is  important  to  remember 
that  the  load  spectra  are  also  normalized  with  respect  to  the  maximum  spectrum  stress 
(usually  maximum  stress=1.0).  Therefore,  it  may  be  necessary  to  adjust  the  Stress 
Multiplication  Factor  (SMF)  to  account  for  the  additional  normalization  of  the  weight 
function  stress  distribution.  The  goal  is  to  ensure  that  the  proper  stresses  are  being 
applied  to  the  crack  face. 

To  enter  a  stress  distribution;  first  go  to  the  Insert  menu;  click  Model,  and  then  select 
Geometry  &  Parameters.  When  the  Model  Geometry  &  Parameters  Dialog  Box 
appears,  cliek  on  the  Stress  Distribution  button.  This  action  calls  the  Weight  Function 
Stress  Distribution  Dialog  Box  (Figure  28),  which  contains  several  control  groups. 

A  user  may  enter  stress  distributions  along  both  ‘a’  and  ‘c’  crack  planes  through  a  table 
of  editable  text  fields  in  the  scrolled  window.  The  default  assumes  that  only  two  sets  of 
data  are  used.  If  the  user  wishes  to  use  more  than  two  sets,  the  desired  number  of  the  sets 
in  the  ‘  Y’  and  ‘X’  directions  must  be  entered.  In  cases  where  the  input  stress  distribution 
does  not  cover  the  extent  of  the  specimen  dimensions,  the  last  input  stress  value  is 
assumed  for  the  remaining  model’s  geometry.  When  a  1-D  problem  is  active,  the  Sy 
part  of  the  stress  distribution  table  is  “grayed  out.” 

At  the  bottom  of  the  dialog,  a  Preview  Window  has  been  placed  in  order  for  a  user  to 
visualize  the  entered,  tabular  stress  distribution.  In  order  to  plot  the  inputted  stress 
distribution  press  the  Apply  button. 

The  two  toggle  buttons  near  the  left  side  of  the  Preview  Window  are  used  to  switch 
between  plots  of  the  ‘S,,’  and  ‘Sy’  stress  distributions. 
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Figure  28.  Weight  Function  Stress  Distribution  Dialog 

The  functionality  of  other  controls  in  Weight  Function  Stress  Distribution  Dialog  Box 
are: 

1 .  NORMALIZE  button:  Normalize  stress  distribution  with  respect  to  the  stress  value 
in  the  crack  origin  row. 

2.  APPLY  button:  Apply  the  current  stress  distribution  to  the  chosen  model. 
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3.  READ  button:  Read  files  with  a  “.wf  s”  extension  ,  which  contain  stress  distribution 
information. 

4.  SAVE  button:  Save  the  current  stress  distribution  to  a  file. 

5.  DISMISS  button:  Dismiss  the  dialog  box. 
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5. 


MODELING  OF  FATIGUE  GROWTH  OF  MULTIPLE  CRACKS 


There  are  several  ways  of  modeling  the  fatigue  growth  of  multiple  cracks.  The  two 
methods  described  briefly  below  are  concerned  with  multiple  cracks  growing  in  the  same 
plane. 

Method  1: 

The  simplest  method  is  to  neglect  the  interference  between  cracks  and  to  model  only  the 
coalescence  of  cracks  as  they  approach  each  other.  The  coalescence  of  cracks  is  often 
model  assuming  the  resultant  dimensions  of  the  final  crack  as  a  sum  of  all  cracks  joining 
each  other  as  shown  in  Figure  29. 


Figure  29.  Coalescence  of  Two  Through-Cracks  in  a  Finite  Width  Plate 


It  is  assumed  in  that  type  of  fatigue  crack  growth  modeling  that  cracks  grow  as  they  were 
single  cracks  with  stress  intensities  at  each  tip  calculated  as  for  a  single  crack  (stage  1). 
They  subsequently  join  each  other  as  soon  as  their  tips  come  to  contact  with  each  other 
and  fi-om  that  point  on  the  two  cracks  are  considered  to  be  transformed  into  one  crack, 
with  the  dimension  being  the  sum  of  appropriate  lengths  of  all  cracks  in  contact  (stage  2). 
In  the  case  of  semielliptical  surface  cracks  it  is  assmned  that  two  or  more  cracks  form  one 
resultant  crack  as  soon  as  their  surface  tips  come  into  contact  with  each  other  (Figure  30). 
The  resultant  length  is  the  sum  of  the  surface  lengths  of  all  cracks  in  contact.  The  depth  is 
assumed  to  be  equal  to  the  depth  of  the  deepest  crack  among  all  cracks  forming  at  a  given 
moment,  the  single  resultant  crack. 
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Figure  30.  Coalescence  of  Two  Semielliptical  Surface  Cracks  in  a  Finite  Thickness  Plate. 
Method  2: 

The  second  method  is  based  on  the  observation  that  fatigue  cracks  join  each  other  very 
fast  as  soon  as  their  cyclic  plastic  zones  at  the  adjacent  tips  come  into  contact.  It  is 
assumed  that  the  cracks  grow  independently  of  each  other  but  they  join  immediately 
when  the  notch  tip  plastic  zones  come  in  contact  as  shown  in  Figure  3 1 . 
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Figure  31.  Coalescence  of  Fatigue  Cracks  Due  to  Sudden  Jump  Over  the  Cyclic  Plastic  Zones 

Coming  in  Contact  with  Each  Other 


Because  the  plastic  zones  contain  highly  damaged  material  it  is  assumed  that  the  crack 
tips  suddenly  jump  over  the  plastic  zones  as  soon  as  they  touch  each  other.  Method  #2  is 
more  conservative  because  it  assumes  that  the  fatigue  crack  growth  over  the  touching 
plastic  zones  occurs  in  one  cycle.  The  resultant  crack  size  is  determined  as  the  sum  of  the 
physical  crack  length  of  each  crack  plus  the  diameter  of  the  cyclic  plastic  zone  of  each 
crack  in  contact.  The  diameter  of  the  cyclic  plastic  zone  in  plane  stress  is  calculated  as: 


J_ 
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The  stress  range  and  the  plastic  zone  size  are  well  defined  in  the  case  of  constant 
amplitude  loading.  However,  in  the  case  of  variable  amplitude  loading  the  stress  range 
have  to  be  defined  by  the  user  depending  on  the  problem  being  analyzed.  It  may  be  the 
largest  range  in  the  loading  history  or  it  can  be  the  average  value  determined  at  given 
crack  size. 


In  the  case  of  the  surface  semielliptical  cracks  it  is  assumed  that  the  crack  join  each  other 
as  soon  as  their  plastic  cones  on  the  surface  come  into  contact.  The  resultant  length  on  the 
surface  is  the  sum  of  the  length  of  each  crack  in  contact  plus  the  diameter  of  each  cyclic 
plastic  zone  coming  into  contact  with  each  other  (Figure  32). 
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The  cracks  are  considered  to  be  “unaware”  of  each  other  until  their  plastic  cyclic  zones 
come  into  contact  with  each  other.  However,  in  order  to  model  the  coalescence  of  fatigue 
cracks,  not  only  their  growth  has  to  be  monitored  continuously  but  their  positions  with 
respect  to  each  other  as  well.  Therefore,  modeling  of  the  growth  of  multiple  fatigue 
cracks  requires  simultaneous  monitoring  of  the  growth  and  position  of  several  cracks  with 
respect  to  a  chosen  global  system  of  coordinates.  The  task  is  intellectually  not  difficult 
but  it  may  pose  a  labor  intensive  programming  task  and  lengthy  numerical  calculations. 
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Figure  32.  Coalescence  of  Two  Semielliptical  Cracks  Over  the  Contacting  Cyclic  Plastic  Zones 
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APPENDIX  A 


73 


NOTATION  FOR  WEIGHT  FUNCTIONS 


A  -  crack  length  for  one-dimensional  cracks  or  crack  depth  (minor  semi¬ 
axis)  for  semielliptical  cracks 

Aj,  Bi  -  coefficients  of  the  linearized  stress  function  in  the  subinterval  “i” 

c  -  crack  length  (major  semiaxis)  for  semielliptical  crack 

Cy  -  coefficients  of  the  integrated  weight  function  mA(x,a)  in  association  with 
the  linearized  stress  field  s(x) 

Dy  -  coefficients  of  the  integrated  weight  function  mB(x,a)in  association  with  the 
linearized  stress  field  s(x) 

K  -  stress  intensity  factor 

Kj  -  mode  I  stress  intensity  factor  (general) 

-  mode  I  reference  stress  intensity  factor  for  the  deepest  point  A  of  a  crack 
under  the  uniform  stress  field 

By,®  -  mode  I  reference  stress  intensity  factor  for  the  surface  point  B  of  a  crack 
xmderthe  uniform  stress  field 

-  mode  I  stress  intensity  factor  for  the  deepest  point  A  on  the  crack  front 
K®  -  mode  I  reference  stress  intensity  factor  for  the  surface  point  B  on  the  crack 

front 

M;  -  parameters  of  the  weight  function  (i  ==  1, 2,  3) 

Mja  -  coefficients  of  the  weight  functions  for  the  deepest  point  A  (i=  1 ,2,3) 

Mjb  -  coefficients  of  the  weight  functions  for  the  surface  point  B  (i=  1 ,2,3) 

m(x,a)  -  weight  function 

mA(x,a)  -  weight  function  for  the  deepest  point  A  of  a  crack 
mB(x,a)  -  weight  function  for  the  surface  point  B  of  a  crack 
S  -  area  under  a  monotonic  curve  m(x,a) 

Sj*  -  area  under  the  linearized  stress  function  s(x)  corresponding  to  the 
subinterval  “i” 

Sj  -  area  under  the  weight  function  curve  m(x,a)  corresponding  to  the  sub- 
interval  “i” 

Q  -  the  semielliptical  crack  shape  factor:  Q  =  1  +1 .464(a/c)*’®^  for  a/c  <  1 
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t  -  wall  thickness  of  a  cylinder  or  plate  thickness 

X  -  the  local,  through  the  thickness  coordinate  along  the  crack  depth 

X  -  coordinate  x  of  the  eentroid  of  the  area  S 

X;  -  coordinate  x  of  the  centroid  of  the  area  Sj  corresponding  the  sub¬ 

interval  “i” 

Xi’  -  coordinate  x  of  the  centroid  of  the  area  S*  corresponding  to  the  sub¬ 
interval  “i” 

a(x)  -  stress  distribution 
a(Xi)  -  value  of  stress  function  at  x  =  Xj 

a(X)  -  value  of  stress  function  at  x  =  X 

ap(x)  -  a  reference  stress  distribution 

-  the  geometric  stress  intensity  correction  factor  for  the  deepest  point  A 

-  the  geometric  stress  intensity  correction  factor  for  the  surface  point  B 
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WEIGHT  FUNCTIONS  FOR  CRACKS  IN  PLATES  &  DISKS 


The  weight  functions  are  given  in  the  form  of  expressions  describing  the  M;  parameters  as 
functions  of  crack  dimensions  and  geometry  of  the  cracked  body.  Given  are  the  range  of 
application,  the  accuracy  and  the  source  of  the  reference  stress  intensity  factors  for  each 
set  of  parameters  M;  including  the  generic  geometry  of  the  crack  body. 


Through-Crack  in  an  Infinite  Plate  Subjected  to  Symmetrical  Loading. 
WEIGHT  FUNCTION  (Figure  33a): 


ml 


(x,a): 


i/27i(a-x) 


1  +  Mj 


1- 


X 


+  M- 


a; 


/  xV 

1-- 

V  aj 


+  M. 


1-^ 


3 


PARAMETERS: 
M,=  0.0698747 
M2=  -0.0904839 
M3=  0.427203 


ACCURACY:  Max.  error  less  than  1%  when  compared  to  exact  solution 

REFERENCES: 

-  weight  function: 

Glinka  G.,  Shen  G..  1991,  “Universal  Features  of  Weight  Functions  for  Cracks  in  Mode 
I,”  Engng.  Fract.  Mech,  Vol.  40,  No.  6,  pp.  1135-1 146. 

-  reference  data: 

Tada  H..  Paris  P.  C.,  Irwin  G.R..  1985,  The  Stress  Analysis  of  Cracks  Handbook,  2nd  ed., 
Paris  Productions  Inc.,  St.  Louis,  Missouri,  pp.  5.1 1-5.1  la. 
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Figure  33.  a)  Through-Crack  in  an  Infinite  Plate  Subjected  to  Symmetrical  Loading;  b)  Central 
Through-Crack  in  a  Finite  Width  Plate  Subjected  to  Symmetrical  Loading 
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Central  Through-Crack  in  a  Finite  Width  Plate  Subjected  to  Symmetric  Loading 


WEIGHT  FUNCTION  (Figure  33b): 


m(x,a)  = 
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RANGE  OF  APPLICATION:  0  <  a/w  <  0.9. 

ACCURACY:  Better  than  1%  when  compared  to  the  reference  solution. 
REFERENCES: 

-  weight function: 

Moftakhar  A..  Glinka  G..  1992,  “Calculation  of  Stress  Intensity  Factors  by  Efficient 
Integration  of  Weight  Functions,”  Engng.  Fract.  Meek,  Vol.  43,  No.  5,  pp.  749-756. 

Glinka  G..  Shen  G.,  1991,  “Universal  Features  of  Weight  Functions  for  Cracks  in  Mode 
I,”  Engng.  Fract.  Meek,  Vol.  40,  No.  6,  pp.  1135-1146. 

-  reference  data: 

Tada  H..  Paris  P.  C..  Irwin  G.R..  1985,  The  Stress  Analysis  of  Cracks  Handbook,  2nd  ed., 
Paris  Productions  Inc.,  St.  Louis,  Missouri,  pp.  2. 1-2.2,  2.34. 
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Sih  G.,  1973,  Handbook  of  Stress  Intensity  Factors,  Institute  of  Fracture  and  Solid 
Mechanics,  Lehigh  University,  Bethlehem,  P.A. 
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Edge  Crack  in  a  Semi-Infinite  Plate 


WEIGHT  FUNCTION  (Figure  34a): 


m(x,a)  = 


^2n{a.-x) 
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+m/ 


1  +  M, 


X 

ay 


1-- 

ay 


3 

^2 


PARAMETERS: 

M,=  0.0719768 
M2=  0.246984 
M3=  0.514465 

ACCURACY:  Max.  error  less  than  1%  when  compared  to  the  reference  solution 
REFERENCES: 

-  weight  function: 

Glinka  G..  Shen  G.,  1991,  “Universal  Features  of  Weight  Functions  for  Cracks  in  Mode 
I,”  Engng.  Fract.  Mech,  Vol.  40,  No.  6,  pp.  1135-1 146. 

-  reference  data: 

Tada  H..  Paris  P.  C..  Irwin  G.R..  1985,  The  Stress  Analysis  of  Cracks  Handbook,  2nd  ed., 
Paris  Productions  Inc.,  St.  Louis,  Missouri,  pp.  8.3-8.3a. 

Sih  G..  1973,  Handbook  of  Stress  Intensity  Factors,  Institute  of  Fracture  and  Solid 
Mechanics,  Lehigh  University,  Bethlehem,  P.A. 
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Edge  Crack  in  a  Finite  Width  Plate 
WEIGHT  FUNCTION  (Figure  34b): 
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RANGE  OF  APPLICATION:  0  <  a/w  <  0.9. 

ACCURACY:  Better  than  1%  when  compared  to  the  reference  solution. 

REFERENCES: 


-  weight function: 

Moftakhar  A..  Glinka  G..  1992,  “Calculation  of  Stress  Intensity  Factors  by  Efficient 
Integration  of  Weight  Fimctions,”  Engng.  Fract.  Meek,  Vol.  43,  No.  5,  pp.  749-756. 

Glinka  G..  Shen  G..  1991,  “Universal  Features  of  Weight  Functions  for  Cracks  in  Mode 
I,”  Engng.  Fract.  Meek,  Vol.  40,  No.  6,  pp.  1135-1 146. 

-  reference  data: 

Tada  H..  Paris  P.  C.,  Irwin  G.R..  1985,  The  Stress  Analysis  of  Cracks  Handbook,  2nd  ed., 
Paris  Productions  Inc.,  St.  Louis,  Missouri,  p.  2.27. 
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Kava  A.  C...  Erdogan  F..  1980,  “Stress  Intensity  Factors  and  COD  in  an  Orthotropic 
Strip,” J  Fracture,  Vol.  16,  pp.  171-190. 


Two  Symmetrical  Edge  Cracks  in  a  Finite  Width  Plate  Subjected  to  Symmetrical 

Loading 


WEIGHT  FUNCTION  (Figure  34c): 

2 


m 


(x,a)  = 
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RANGE  OF  APPLICATION:  0  <  a/w  <  0.9. 

ACCURACY:  Better  than  1%  when  compared  to  the  reference  solution. 

REFERENCES: 

-  weight function: 

Moftakhar  A..  Glinka  G..  1992,  “Calculation  of  Stress  Intensity  Factors  by  Efficient 
Integration  of  Weight  Functions, “  Engng.  Fract.  Meek,  Vol.  43,  No.  5,  pp.  749-756. 

Glinka  G..  Shen  G.,  1991,  “Universal  Features  of  Weight  Functions  for  Cracks  in  Mode 
I,”  Engng.  Fract.  Meek,  Vol.  40,  No.  6,  pp.  1135-1146. 

-  reference  data: 

Tada  H..  Paris  P.  C..  Irwin  G.R..  1985,  The  Stress  Analysis  of  Cracks  Handbook,  2nd  ed., 
Paris  Productions  Inc.,  St.  Louis,  Missouri,  p.  2.31. 
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Shallow  Semielliptical  Surface  Crack  in  a  Finite  Thickness  Plate  (a/c  <  1) 


WEIGHT  FUNCTION  (Figure  35)  -  Deepest  Point  A: 
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RANGE  OF  APPLICATION:  0  =  a/t  =0.8  and  0  =  a/c  =  1.0 
ACCURACY:  Better  than  3%  when  compared  to  the  FEM  data. 

REFERENCES: 

-  weight  function: 

Wang  X..  Lambert  S.  B..  1995,  “Stress  Intensity  Factors  for  Low  Aspect  Ratio 
Semielliptical  Surface  Cracks  in  Finite-Thickness  Plates  Subjected  to  Non-uniform 
Stresses,”  Engng.  Fract.  Mech,  Vol.  51,  No.  4,  pp.  517-532. 

Shen  G..  Glinka  G..  1991,  “Weight  Functions  for  a  Surface  Semi-Elliptical  Crack  in  a 
Finite  Thickness  Plate,”  Theor.  Appl.  Fract.  Mech,  Vol.  15,  No.  2,  pp.  247-255  (this 
paper  contains  the  older  version  of  the  weight  function). 

Shen  G..  Plumtree  A..  Glinka  G..  1991,  “Weight  Function  for  the  Surface  Point  of  Semi- 
Elliptical  Surface  Crack  in  a  Finite  Thickness  Plate,”  Engng.  Fract.  Mech,  Vol.  40,  No. 
1,  pp.  167-176  (this  paper  contains  the  older  version  of  the  weight  function). 
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-  reference  data: 

Shiratori  M.,  Mivoshi  T..  Tanikawa  K.,  1986,  “Analysis  of  Stress  Intensity  Factors  for 
Surface  Cracks  Subjected  to  Arbitrarily  Distributed  Surface  Stress  (2nd  Report,  Analysis 
and  Application  of  Influence  Coefficients  for  Flat  Plates  with  a  Semielliptical  Surface 
Crack),“  Trans.  JSME,  Vol.  52,  pp.  390-398. 

Gross  B..  Srawlev  J.  E..  1965,  “Stress  Intensity  Factor  for  Single-Edge-Notch  Specimens 
in  Bending  or  Combined  Bending  and  Tension  by  Boundary  Collocation  of  a  Stress 
Function,”  NASA  TN,  D-2603. 

Gross  B..  Srawlev  J.  E..  1 966,  Plane  Strain  Crack  Toughness  Testing  of  High  Strength 
Metallic  Materials,  ASTMSTP  410. 

Raiu  I.  S.,  Newman  J.  C..  1979,  “Stress  Intensity  Factors  for  a  Wide  Range  of 
Semielliptical  Surface  Cracks  in  Finite  Thickness  Plate,”  Engng.  Fracture  Meek,  Vol. 
11,  pp.  817-829. 

Newman  J.  C..  Raiu  I.  S..  1981,  “An  Empirical  Stress  Intensity  Factor  Equation  for  the 
Surface  Crack,”  Engng.  Fracture  Meek,  Vol.  15,  pp.  185-192. 


Figure  35.  Shallow  Semielliptical  Surface  Crack  in  a  Finite  Thickness  Plate 
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Deep  Semielliptical  Surface  Crack  in  a  Finite  Thickness  Plate  (a/c  >  1) 


WEIGHT  FUNCTION  (Figure  36)  -  Deepest  Point  A: 
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WEIGHT  FUNCTION  (Figure  36)  -  Surface  Point  B 
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f-1 

UJ 

and 


F,= 


Do  +D, 


U, 


+  D. 


D„  =1.11855-0.2065 


vt; 

+  0.078171 


VC7 


D,  =1.15312-0.98743 


>Y 

^cj 


+  0.23315 


D2  =0.2246-0.47841 


-1  +  0.18641 
VC 


a 

-c; 


a 

VC 


RANGE  OF  APPLICATION:  0  =  a/t  =  0.8  and  0.6  =  a\c  =  2.0 
ACCURACY:  Better  than  2%  when  compared  to  the  reference  FEM  data. 
REFERENCES: 


-  weight  function: 

Wang  X..  1995,  “Weight  Functions  for  High  Aspect  Ratio  Semielliptical  Surface  Cracks 
in  Finite-Thickness  Plates,”  University  of  Waterloo,  to  be  published. 

-  reference  data: 

Shiratori  M..  Mivoshi  T..  Tanikawa  K..  1986,  “Analysis  of  Stress  Intensity  Factors  for 
Surface  Cracks  Subjected  to  Arbitrarily  Distributed  Surface  Stress  (2nd  Report,  Analysis 
and  Application  of  Influence  Coefficients  for  Flat  Plates  with  a  Semielliptical  Surface 
Crack),”  Trans.  JSME,  Vol.  52,  pp.  390-398. 
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Quarter-Elliptical  Corner  Crack  in  a  Finite  Thickness  Plate 


WEIGHT  FUNCTION  (Fig.  37)  -  Deepest  Point  A: 


PARAMETERS: 

M«=3 

M„=^(2Y,-Y.)4 

where: 


Q  =  1  +  1.464 


ra^  raV  .  raV  A  faV 

y„=Ao+a(-J+a,[-J  +A3[-J  .a.^-J 

fa')  TaV  faV 

A,  =1.041  +  0.016  -  +0.186  -  -0.111  - 

“  U;  ycj 

fa^  fa']^  /^aV 

A,  =-0.599  +  1.953  - 1.310  -  -0.028  - 

'  ycj  \cj 

fa'\  TaV  TaV 

A,  =4.972-13.216  -  +6.747  -  +1.918  - 

^  Vc7  vcj  Vcj 

fa\  TaV  TaV 

A,  =-1.293  +  1.857  -  +12.906  -  -13.441  - 

^  {cj  VcJ  VcJ 

A,  =-0.572  +  3.073f-Vl0-797f-l  +8.393f-l 
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Y,  =B„  +  B,  -  +BJ-  +I 


B„  =0.500-0.323  -  +0.213 


B,  =-0.507  +  1.373  -  -0.74' 


Bj  =3.468-9.028  -  +6.349 


B3  =-1.359  +  1.731  -  +5.35 

U  J 


B.  =-0.162  +  2.977  -  -8.25 


WEIGHT  FUNCTION  (Figure  37)  -  Surface  Point  B 


mg 


(x,a) 


1  +  M 


IB 


f-1 

UJ 


1 


+  M 


2B 


'-T 


+  M33 


f-1 

UJ 


3 

2 


PARAMETERS: 


M„=^(5F,-3F,)-8 

M„=^(2F.-3F,)+15 

M,,=^(lOF,-7F,)-8 


where: 


Q  =  l  +  1.464 


/  \1.65 

^  a  ^ 


Fo  = 


Co+C, 


vCy 

fi 


+  C 


/^o^ 


{tj 


+  C, 


(-] 

KiJ 


u. 


Co  =3.340-4.495 

ycj 

C,  =0.2318-0.2261 


'^a'' 

VC; 


vc; 


C2  =22.95-100.9 


Vc; 


C3  =-39.16  +  194.1 


vc; 


C4  =30.80-142.91 


vc; 


+  3.016| 
-1.6581 
+  152.2| 
-302.01 
+  212.6! 


-0.7278  - 

vcy 

2 

+  1.5041 


Vcj 


UY 

Vc> 


-72.92 


vc; 


UT 

VC; 
^a^- 


vc; 

+  145.9f-^ 

VcJ 


vcy 


-99.92 


cj 


and 
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Fi  = 


Dq  +D, 


V  Ly 


+  D. 


+  D. 


Do  =2.831-3.8401 


^a^ 


vty 
+  2.4771 


+  D/ 


'-T 

VCy 


(-] 

'oV 


^a^ 


vc; 


D.  =4.600-20.498 


D2  =-4.019  +  15.057 


+  29.001 

UJ 

UJ 

-0.511  - 

ycj 

2 

-13.2261 


-12.624 


^a^' 


D3  =9.682-15.932 


a 

VC 


vc; 

-8.8481 


vc; 


+  2.6771 

+  I3.910I 


vc; 


Vc; 


D.  =-1.141-9.176 


Vcy 


vc; 

+  30.2281 


^a^^ 


vc; 


•19.195  - 


vc; 


cj 


RANGE  OF  APPLICATION:  0  =  a\t  =0.8  and  0.2=  a\c  =  1.0 
ACCURACY:  Better  than  1.5%  when  compared  to  the  FEM  data. 
REFERENCES: 


-  weight function: 

Zheng  X.  J..  Glinka  G..  Dubev  R.  N..  1994,  “Stress  Intensity  Factors  and  Weight 
Functions  for  a  Comer  Crack  in  a  Finite  Thickness  Plate,”  Engng.  Fracture  Mech.,{io  be 
published,  1996). 

-  reference  data: 

Shiratori  M..  Mivoshi  T..  1985,  “Analysis  of  Stress  Intensity  Factors  for  Surface  Cracks 
Subjected  to  Arbitrarily  Distributed  Surface  Stress  (Analysis  and  Application  of  Data 
Base  of  Influence  Coefficients  Kjj),”  Proc.  Third  Conf.  Fracl.  Meek,  Soc.  Materials 
Science  Japan,  pp.  82-86. 

Murakami  Y..  et  ak.  1992,  Stress  Intensity  Factors  Handbook,  Vol.  3,  Pergamon  Press, 
pp.  591-597. 

Raiu  L  S.,  Newman  J.  C..  1988,  “Stress  Intensity  Factors  for  Comer  Cracks  in 
Rectangular  Bars,“  in  Fracture  Mechanics:  Nineteenth  Symposium,  ASTM  STP  969, 
T.A.Cmse,  ed.,  pp.  43-55. 
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Figure  37.  Quarter-Elliptical  Corner  Crack  in  a  Finite  Thickness  Plate 
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Radial  Edge  Crack  in  a  Circular  Disk 
WEIGHT  FUNCTION  (Figure  38): 


PARAMETERS: 


f  a  ^  /"aV 

-0.04732  +  0.49586  - 1.94141  —  +3.96175  - 

vdJ 

/  \4  / 

-4.56109  —  +2.83026  —  -0.74080  — 

VDJ  VDy  VDJ 

f  f  a^^  f  a\^ 

-0.58602  +  8.48276  —  -30.90993  —  +112.21153  — 

vdJ  vdJ  I^dJ 

/  \4  /  n5  f'  f 

-280.25303  —  +428.48183  —  -356.66155  —  +125.34267  — 

VDJ  VDJ  i^Dj 

f  a\  f  a\^  f  a\^ 

-1.09836  +  3.06605  —  +16.85709  —  -48.14897  — 

vdJ  vdJ  vdJ 

/  \4  f  f  f 

+  54.61627  -  +6.91042  -  -61.35817  —  +36.1270  — 

vdJ  vdJ  vDj  vdJ 

RANGE  OF  APPLICATION:  0  <  a/D<  0.9. 

ACCURACY:  Better  than  1.5%  compared  to  the  reference  solution. 

REFERENCES: 

-  weight  function: 

Kiciak  A..  1995,  University  of  Waterloo,  to  be  published. 

Glinka  G..  Shen  G.,  1991,  “Universal  Features  of  Weight  Functions  for  Cracks  in  Mode 
I,”  Engng.  Fract.  Meek,  Vol.  40,  No.  6,  pp.  1 135-1 146  -  this  paper  contains  the  original 
version  of  the  weight  function  with  coefficients  Mj  given  in  tabular  form. 

-  reference  data: 

Wu  X.  R..  Carlsson  J.  A..  1991,  Weight  Functions  and  Stress  Intensity  Factor  Solutions, 
Pergamon  Press. 
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WEIGHT  FUNCTIONS  FOR  CRACKS  IN  CYLINDERS 


The  weight  functions  are  given  in  the  form  of  expressions  describing  the  M;  parameters  as 
functions  of  crack  dimensions  and  geometry  of  the  cracked  body.  Given  are  the  range  of 
application,  the  accuracy  and  the  source  of  the  reference  stress  intensity  factors  for  each 
set  of  parameters  Mj  including  the  generic  geometry  of  the  crack  body.  The  geometry  of  a 
cylinder  is  characterized  by  the  ratio  of  the  external  to  internal  radius,  Ro/Rf. 
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Internal  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder  (R„/Ri=2.0) 


WEIGHT  FUNCTION  (Figure  39)  -  Deepest  Point  A: 


mA(x,a)  = 


(  v\2  {  vV  f  XV 

1  +  M,a  1--  +M3A  1—  +M3A  1  — 

vaj  vay  lay 


PARAMETERS: 


Mm  =3 


where: 


Q  =  1  +  1.464 


and  t  =  R..  -  R, 


*  ,  faV  A  faV 

Yo=Ao+Ai  -  +A,  -  +A3  - 

Uy  vty  \ij 

A„=1.12-0.207f-Vo.l53f-T+1.305f-l  -2.007f-l  +0.933f- 


A,  =-0.111-7.205f-V36.455f-l  -83.649f-l  +89.74lf-l  -35.219f-l 

'  UJ  vc^  Uy  vcy  vcy 

Aj  =1.498  +  20.265  ^-j-132.935|^^j  +323.535|^^  -343.92oj^^j  +131.532^-^ 

A3  =-0.140  +  18.828f-V85.243f-l  +118.94lf-l  -52.084f-] 

^  Ui  UJ  Uy  \^J 
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2 


3 


4 


WEIGHT  FUNCTION  (Figure  39)  -  Surface  Point  B: 

1  A/T  f  IV/r  f  A/f 

1  +  M,b  —  +M2B  —  +M 

U;  U; 

PARAMETERS 
M„=^(2F.-5F,)-8 

1  Sir 

M3b='^(3F,-10F,)-8 


Q  =  1  +  1.4641 


/  \l-65 

'  a  1 


and  t  =  R„-R: 


Fo  = 


Cp  +C, 


^-1 

ty 


+  C. 


+  C.. 


viy 


a  =5.923-20.55 


Cj  =-3.607  +  11.6861 


TaV 

+  36.937  - 


vc; 

-31.6341 


f-T 

U; 


Cj  =19.14-72.902 


C3  =9.586-64.389 


v'-y 

^-1 

Vcy 

^-1 

Vcy 


14.138 
+  112.643] 
+  151.4491 


^a^' 


^-T 

UJ 

^-T 


+  9.935 

-86.9041 


+  10.37 
-3.7741 


^a^^ 


vc; 

'  a ' 


Vcj 

a 


vc; 

+  28.149 


-144.822 


a 

vc; 


(-] 

3 

+  48.124 

f-1 

UJ 

UJ 

and 


4 

Do+D,  - 

+  D2 

+  F>3  — 

— 

UJ 

U; 

vcy 

F,= 


Dn  =0.687-1.821 


Uy 

y 


+  2.718 


^-1 

Vcy 


-1.981 


+  0.567 


v'^y 


a^ 

.cy 


Di  =-1.797  +  9.399 


-  -19.1951  - 


D2  =  8.504-42.6081 


vc; 

^-1 


+  87.828 


Kcj 

^-T 

'.cy 


+  17.88lf-^ 

vcy 


■82.735 


VoJ 


-6.244 

Vcy 

fa^ 
+  29.055  - 
vc; 


D,  =  -0.666  +  6.54 


-21.603 


Vcy 


y< 


+  27.036 


'.c 


-  -11.326  - 


RANGE  OF  APPLICATION:  R^/R;  =  2.0,  0  =  a/t  =  0.8,  a/c  =  0  and  0.2  =  a/c  =  1.0 
ACCURACY:  Better  than  3%  when  compared  to  the  FEM  data. 

REFERENCES: 


-  weight  function: 

Zheng  X.  J..  Glinka  G..  1995,  “Weight  Functions  and  Stress  Intensity  Factors  for 
Longitudinal  Semi-Elliptical  Cracks  in  Thick-Wall  Cylinders,”  J.  Press.  Vess.Technol., 
ASME,  vol.  117, 1995,  pp.  383  -  389 

-  reference  data: 

Mettu  S.  R..  Raiu  I.  S..  Forman  R.  G..  1988,  “Stress  Intensity  Factors  for  Part-Through 
Surface  Cracks  in  Hollow  Cylinders,”  NASA  Technical  Report,  No.  JSC  25685  LESC 
30124. 
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Andrasic  C.  P..  Parker  A.  P..  1984,  “Dimensionless  Stress  Intensity  Factors  for  Cracked 
Thick  Cylinder  Under  Polynomial  Crack  Face  Loading,”  Engng.  Fract.  Meek,  Vol.  19, 
1984,  pp.  187-193. 

Murakami  Y..  et  al..  1986,  Stress  Intensity  Factor  Handbook,  Vol.  2,  Pergamon  Press,  pp. 
309-317. 


Figure  39.  Internal  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder 
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Internal  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder  (Rj,/Rj-1.5) 


WEIGHT  FUNCTION  (Figure  39)  -  Deepest  Point  A: 


mA(x,a)=  ,  .  -T 

^27t(a-x) 


1  +  M 


lA 


^  xVa 


V  a; 


+  M 


2A 


1-- 

V  a; 


+  M 


3A 


1- 
V  a 


xP 


PARAMETERS: 


Mm  =3 


where: 


Q  =  1  +  1.464 


(-) 

UJ 

a 


1.65 


and  t  =  -  Rj 


+  A- 


Yq  =  Aq  +  A] 

Aq  =1.044 +  0.07  exp 
A,  =  -0.433  +  0.665  exp 
Aj  =0.711 +  1.161  exp 
Aj  =-0.179 +  1.46  exp 


+  A: 


a. 


5.051 


^  a. 

VC, 


/a' 

-3.393  - 
VC, 

-3.386f-^ 

Uy 


-4.165 


^a' 

VC, 


and 
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N»  I  OJ 


Q  =  1  +  1.464 


1.65 


Fo  = 


Cq  +C] 


yoj 


and  t  =  Rq  -  Rj 


+  C, 


t 


a 

vt; 


a^ 

.c; 


Co  =  0.972  +  5. 163  exp 


-5.061 


1  +  1.568 
yo) 


Cl  = -0.199 -10.239exp| 
Cj  =0.119  +  8.784exp 
C3  = -0.104 +  28. 133exp| 


-46.053 


\.oJ 


-4.0091 


-4.081 


-1  +  1.092 
KcJ 


-9.959f- 

U; 


-9.817  - 
vcy 


and 


F,= 


Do  +Di 


(-) 


+  D. 


Dn  =1.033-4.842 


Vt. 

+  9.7081 


f-T 


-8.397 


D,  =-3.448  +  24.231 


Kcj 

•50.221 


(-] 


+  2.690 


'  a ' 


D,  =6.535 -30.622^^^ 


\ 


D3  =2.243-21.6771  - 


a 

cj 


+  45.644] 
+  65.5461 


VC7 

\cJ 

a 

YcJ 


+  42.498 

ycj 

I  -25.05(iJ 


vc; 

-13.099 

/ 


a 

VC. 


N‘’ 

- 

vc; 

4 


7oV 


-76.555 


a 

vCy 


+  3.636 

+  30.4331 


vc; 


RANGE  OF  APPLICATION:  Rg/Ri  =  1.5,  0  =  a/t  =  0.8,  a/c  =  0  and  0.2  =  a/c  =  1.0 
ACCURACY:  Better  than  3%  when  compared  to  the  FEM  data. 

REFERENCES: 


-  weight  function: 

Zheng  X.  J..  Kiciak  A..  Glinka  G..  1996,  “Determination  of  Weight  Fimctions  and  Stress 
Intensity  Factors  for  Semielliptical  Cracks  in  Thick-wall  Cylinders,”  Nucl  Engng. 
Design,  to  be  published. 
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-  reference  data: 

Mettu  S.  R..  Raiu  I.  S..  Fonnan  R.  G..  1988,  “Stress  Intensity  Factors  for  Part-Through 
Surface  Cracks  in  Hollow  Cylinders,”  NASA  Technical  Report,  No.  JSC  25685  LESC 
30124. 

Andrasic  C.  P.,  Parker  A.  P..  1984,  “Dimensionless  Stress  Intensity  Factors  for  Cracked 
Thick  Cylinder  Under  Polynomial  Crack  Face  Loading,”  Engng.  Fract.  Meek,  Vol.  19, 
1984,  pp.  187-193. 

Murakami  Y.,  et  al..  1986,  Stress  Intensity  Factor  Handbook,  Vol.  2,  Pergamon  Press,  pp. 
309-317. 
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Internal  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder  (R„/Rj=1.25) 


WEIGHT  FUNCTION  (Figure  39)  -  Deepest  Point  A: 


0  (  {  xV  {  X12 

mA(x,a)  =  -^  --  1  +  M,a  1--  +^2A  1--  +^3a|i-- 

Jlnia-x)  \  \  V  ay 


PARAMETERS: 


M«=3 


where: 


Q  =  1  +  1.464 


and  t  =  R„-R, 


faV  faX 

Yq  =  Aq  +  A]  -  +  A2  -  +A3  - 

vty  \i J  \ij 


Ao  =1.010  + 0.0998exp  -13.15  - 

\c 


A,  =  0.055 +  0.366exp  -31.17  - 

V  c 


A,  =  -0.057  +  3. 269 exp  -  3.859 


A3  =-0.149 +  0.061  exp  1.3541 - 
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2 


3 


1.65 


Q  =  l  +  1.464 


a 


Fo  = 


VC 

f 


and  t  =  Ro-Ri 


Co+C, 


i] 


t; 


+c. 


vty 


+  C, 


Cn  =5.566-19.583 


C,  =-1.75  +  9.514 


Vc7 

^  a^ 


vc; 


+  37.3351 
-16.6181 


'-T 


vc; 


-33.705 


vc; 


vc; 


+  11.507 


^a^^ 


Vc; 


+  10.44 


=12.497-49.067 


vc; 

+  72.591 


f-1 

UJ 


-1.616 


vc; 


^a^^ 


Cj  =3.468-29.49! 
and 


vc; 

+  83.789 


vc; 


-45.216  - 

U; 

.3 


+  9.55 


f-1 

UJ 


-93.289 


vc; 


+  35.507 


vc; 


F,= 


■ 

fa') 

1  fa.') 

2 

faV 

(-] 

Do  +D, 

+  D3 

7 

U  1  1 

Uy 

1  \tJ 

UJ 

Dn  =0.486-0.879 


vc; 


D,  =-0.533  +  2.626  - 
Vcy 


vc; 


D2  =4.116-15.9851 


+  1.161 
-3.412| 
+  22.3581 


-0.793f-l  +0.212f-^ 


vc; 


vc; 

+  O.999I 


vc; 

+  0.333i 


'  a  ' 


VC7 


vcy 
-12.2351 


vc; 


vc; 

+ 1 .826| 


r  N'* 
'  a  ' 


vc; 


D3  =0.569-6.6051 


vc; 


+  21.548 


UJ 

VCy 


-26.37 


Vc7 


+  10.853 


vc; 

r 

’  a 
Vcy 


RANGE  OF  APPLICATION:  Ro/R;  =  1 .25, 0  =  a/t  =  0.8,  a/c  =  0  and  0.2  =  zlc=\S) 
ACCURACY:  Better  than  3%  when  compared  to  the  FEM  data. 

REFERENCES: 


-  weight function: 

Zheng  X.  J.  CTlinka  G..  Dubev  R.  N..  1995,  “Calculation  of  Stress  Intensity  Factors  for 
Semielliptical  Cracks  in  a  Thick-Wall  Cylinder,”  Int.  J.Pres.Ves.  &  Piping,  Vol.  62,  pp. 
249-258. 


Ill 


-  reference  data: 

Mettu  S.  R-  Raiu  I.  S..  Forman  R.  G..  1988,  “Stress  Intensity  Factors  for  Part-Through 
Surface  Cracks  in  Hollow  Cylinders,”  NASA  Technical  Report,  No.  JSC  25685  LESC 
30124. 

Andrasic  C.  P..  Parker  A.  P..  1984,  “Dimensionless  Stress  Intensity  Factors  for  Cracked 
Thick  Cylinder  Under  Polynomial  Crack  Face  Loading,”  Engng.  Fract.  Mech,  Vol.  19, 
1984,  pp.  187-193. 

Murakami  Y..  et  al..  1986,  Stress  Intensity  Factor  Handbook,  Vol.  2,  Pergamon  Press,  pp. 
309-317. 
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Internal  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder  (R^/Rj— 1.1) 
WEIGHT  FUNCTION  (Figure  39)  -  Deepest  Point  A: 


PARAMETERS: 


f 

Q  =  1  +  1.464  -  and  t  =  Ro-Ri 

vc; 

/aV  /^aV 

Yo  =  Ag+Aj  -  +A2  - 

vty  Vty 

f^\  /aV  /^a\ 

Ao  =1.1449-0.6699  -  +1.0464  -  -0.5202  - 

v^y  V^y  v^y 

/^aV 

A,  =3.84-10.531  -  +6.931  - 

vcy  Vcy 


and 
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/aV 

Y,=Bo  +  B.  -  +B,  - 

V  t  i  \i  j 


B„  =0.4732-0.4967  -  +0.7576  -  -0.4417  - 


B,  =2.415-6.901  -  +5.928  -  -1.291  - 


Bj  =-6.251  +  13.282  -  -8.097  -  + - 5^^ 

\.Cy  yCJ  f  “ 

0.090+  - 


WEIGHT  FUNCTION  (Figure  39)  -  Surface  Point  B: 


ma(x,a)  =  ^  1  +  M,a  -  ■+M2B  -  - 


PARAMETERS: 


M,b=^(5F,-3F„)-8 

1 

Mib=^(2Fo-3F,)+15 

M3,=^(10F,-7F„)-8 


where: 
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RANGE  OF  APPLICATION:  VRi  =  1.1,  0  =  a/t  =  0.8,  and  0  =  a/c  =  1.0 
ACCURACY:  Better  than  5%  when  compared  to  the  FEM  data. 

REFERENCES: 

-  weight  function: 

Wang  X.  J.-  Lambert  S.  B..  1996,  “Stress  Intensity  Factors  and  Weight  Functions  for 
Longitudinal  Semielliptical  Surface  Cracks  in  Thin  Pipes,”  Int.  J.Pres.  Ves.  &  Piping, 
Vol.  65,  pp.  75  -87. 

Shen.  G..  Liebster  T.  D..  Glinka  G.  ,1993,  “Calculation  of  Stress  Intensity  Factors  for 
Cracks  in  Pipes,”  Proc.  of  the  12th  Int.  Conf  on  Offshore  Meek  and  Arctic  Engng., 
ASME,  M.  Salama  et  al,  ed.,  Vol.  III-B,  Materials  Engineering,  pp.847-854.  (this  paper 
contains  the  older  version  of  the  weight  function,  limited  to  0  =  a/t  =  0.8  and  0.2  =  a/c  = 
1.0) 
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-  reference  data: 

Raiu  I.  S..  Newman  J.  C..  1982,  “Stress  Intensity  Factors  for  Internal  and  External 
Surface  Cracks  in  Cylindrical  Vessels,”/.  Press.  Ves.  Technol.,  Vol.  104,  pp.  293-298. 

Andrasic  C.  P.,  Parker  A.  P..  1984,  “Dimensionless  Stress  Intensity  Factors  for  Cracked 
Thick  Cylinder  Under  Polynomial  Crack  Face  Loading,”  Engng.  Fract.  Meek,  Vol.  19, 
1984,  pp.  187-193. 


External  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder  (Ro/Rj-l.O) 
WEIGHT  FUNCTION  (Figure  40)  -  Deepest  Point  A: 


m 


■y]2n{a-x) 


f 


1  +  M 


lA 


1 


1- 

V  a; 


+  M 


2A 


1-- 


V 


+  M 


3A 


'■J 


PARAMETERS: 


24 


M,,=3 


where: 


Q  =  1  +  1.464 
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WEIGHT  FUNCTION  (Figure  40)  -  Surface  Point  B: 


(  \  ^  1  ^/r  r  Ayr  A/r  r 

mB(x,a)  =  -=  1  +  M,b  -  +M2B  -  +M3B  - 

^J^tx  \aj  yaj  yaj 


PARAMETERS: 


B,  =0.02433 -0.08062 exp  24.45619 
Bj  =0.12875 +  0.46420  exp  14.79057 
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RANGE  OF  APPLICATION:  R/R^  =  2.0,  0  =  a/t  =  1.0,  a/c  =  0  and  0.2  =  a/c  =  1 .0 
ACCURACY:  Better  than  3%  when  compared  to  the  FEM  data. 

REFERENCES: 
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-  weight  function: 

Kiciak  A..  Glinka  G..  Bums  D.  J..  1996,  “Weight  Functions  and  Stress  Intensity  Factors 
for  Longitudinal  Semielliptical  Cracks  in  Thick-wall  Cylinders,”  J.  Press.  Vess.  Technol. , 
ASME  (to  be  published) 

-  reference  data: 

Mettu  S.  R..  Raiu  I.  S.,  Forman  R.  G..  1988,  “Stress  Intensity  Factors  for  Part-Through 
Surface  Cracks  in  Hollow  Cylinders,”  NASA  Technical  Report,  No.  JSC  25685  LESC 
30124. 

Andrasic  C.  P..  Parker  A.  P..  1984,  “Dimensionless  Stress  Intensity  Factors  for  Cracked 
Thick  Cylinder  Under  Polynomial  Crack  Face  Loading,”  Engng.  Fract.  Meek,  Vol.  19, 
1984,  pp.  187-193. 

Murakami  Y..  et  al..  1986,  Stress  Intensity  Factor  Handbook,  Vol.  2,  Pergamon  Press,  pp. 
309-317. 


Figure  40.  External  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder 


External  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder  (Ro/Rj-l.S) 
WEIGHT  FUNCTION  (Figure  40)  -  Deepest  Point  A: 
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and 
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WEIGHT  FUNCTION  (Figure  40)  -  Surface  Point  B: 
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PARAMETERS: 


M„=;^(2F.-5F,)-8 

Mjb=^(3F,-F.)+15 

M3b=^(3F,-10F,)-8 


where: 
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F]  -  Do  +  D, 
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RANGE  OF  APPLICATION:  R/R^  =  1.5,  0  =  a/t  =  1.0,  a/c  =  0  and  0.2  =  a/c  =  1.0 
ACCURACY:  Better  than  3%  when  compared  to  the  FEM  data. 

REFERENCES: 

-  weight function: 

Kiciak  A..  Glinka  G..  Bums  D.  J..  1996,  University  of  Waterloo,  to  be  published. 

-  reference  data:  ^ 

Mettu  S.  R..  Raiu  I.  S..  Forman  R.  G..  1988,  “Stress  Intensity  Factors  for  Part-Through 
Surface  Cracks  in  Hollow  Cylinders,”  NASA  Technical  Report,  No.  JSC  25685  LESC 
30124. 

Andrasic  C.  P..  Parker  A.  P..  1984,  “Dimensionless  Stress  Intensity  Factors  for  Cracked 
Thick  Cylinder  Under  Polynomial  Crack  Face  Loading,”  Engng.  Fract.  Meek,  Vol.  19, 
1984,  pp.  187-193. 

Murakami  Y..  et  al..  1986,  Stress  Intensity  Factor  Handbook,  Vol.  2,  Pergamon  Press,  pp. 
309-317. 
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External  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder  (Rp/Rj— 1.25) 


WEIGHT  FUNCTION  (Figure  40)  -  Deepest  Point  A: 
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PARAMETERS: 
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RANGE  OF  APPLICATION:  Ro/R;  =  1.25,  0  =  a/t  =  1.0,  a/c  =  0  and  0.2  =  a/c  =  1.0 
ACCURACY:  Better  than  3%  when  compared  to  the  FEM  data. 

REFERENCES: 


-  weight  function: 

Kiciak  A..  Glinka  G..  Bums  D.  J..  1996,  University  of  Waterloo,  to  be  published. 

-  reference  data: 

Mettu  S.  R..  Raiu  I.  S..  Forman  R.  G..  1988,  “Stress  Intensity  Factors  for  Part-Through 
Surface  Cracks  in  Hollow  Cylinders,”  NASA  Technical  Report,  No.  JSC  25685  LESC 
30124. 

Andrasic  C.  P..  Parker  A.  P..  1984,  “Dimensionless  Stress  Intensity  Factors  for  Cracked 
Thick  Cylinder  Under  Polynomial  Crack  Face  Loading,”  Engng.  Fract.  Meek,  Vol.  19, 
1984,  pp.  187-193. 

Murakami  Y.,  et  al..  1986,  Stress  Intensity  Factor  Handbook,  Vol.  2,  Pergamon  Press, 
pp.  309-317. 
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External  Axial  Semielliptical  Surface  Crack  in  a  Thick  Cylinder  (Ro/Ri=l.l) 


WEIGHT  FUNCTION  (Figure  40)  -  Deepest  Point  A: 
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WEIGHT  FUNCTION  (Figure  40)  -  Surface  Point  B: 
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RANGE  OF  APPLICATION:  =  1.1, 0  =  a/t  =  1.0,  and  0  =  a/c  =  1.0 

ACCURACY:  Better  than  5%  when  compared  to  the  FEM  data. 
REFERENCES: 


-  weight  function: 

Wang  X.  J..  Lambert  S.  B..  1996,  “Stress  Intensity  Factors  and  Weight  Functions  for 
Longitudinal  Semielliptical  Surface  Cracks  in  Thin  Pipes,”  Int.  J.Pres.Ves.  &  Piping, 
Vol.  65,  pp  75-87. 

-  reference  data: 

Raiu  I.  S..  Newman  J.  C..  1982,  “Stress  Intensity  Factors  for  Internal  and  External 
Surface  Cracks  in  Cylindrical  Vessels,”  J.  Press.  Ves.  TechnoL,  Vol.  104,  pp.  293-298. 
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Andrasic  C.  P..  Parker  A.  P.,  1984,  “Dimensionless  Stress  Intensity  Factors  for  Cracked 
Thick  Cylinder  Under  Polynomial  Crack  Face  Loading,”  Engng.  Fract.  Meek,  Vol.  19, 
1984,  pp. 187-193. 
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